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ABSTRACT 

We combine stellar metallicity and stellar mass estimates for a large sample of galaxies 
drawn from the Sloan Digital Sky Survey Data Release Two (SDSS DR2) spanning 
wide ranges in physical properties, in order to derive an inventory of the total mass 
of metals and baryons locked up in stars in the local Universe. Physical parameter 
estimates are derived from galaxy spectra with high signal-to-noise (S/N) ratio (of 
at least 20). Coadded spectra of galaxies with similar velocity dispersions, absolute 
r-band magnitudes and 4000A-break values are used for those regions of parameter 
space where individual spectra have lower S/N. We estimate the total density of metals 
pz and of baryons in stars and, from these two quantities, we obtain a mass- and 
volume-averaged stellar metallicity of (Z*) — 1.04±0.14 Zq, i.e. consistent with solar. 
We also study how metals are distributed in galaxies according to different properties, 
such as mass, morphology, mass- and light- weighted age, and we then compare these 
distributions with the corresponding distributions of stellar mass. We find that the bulk 
of metals locked up in stars in the local Universe reside in massive, bulge-dominated 
galaxies, with red colours and high 4000A-break values corresponding to old stellar 
populations. Bulge-dominated and disc-dominated galaxies contribute similar amounts 
to the total stellar mass density, but have different fractional contributions to the mass 
density of metals in stars, in agreement with the mass-mctallicity relation. Bulge- 
dominated galaxies contain roughly 40 percent of the total amount of metals in stars, 
while disc-dominated galaxies less than 25 percent. Finally, at a given galaxy stellar 
mass, we define two characteristic ages as the median of the distributions of mass 
and metals as a function of age. These characteristic ages decrease progressively from 
high-mass to low-mass galaxies, consistent with the high formation epochs of stars in 
massive galaxies. 
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1 INTRODUCTION 

Constraining the star formation and chemical evolution his- 
tories of galaxies is one of the fundamental goals in observa- 
tional cosmology. The evolution of the global star formation 
rate (SFR) has to map the evolution over cosmic time of 
its products, i.e. the baryonic and metal content of the Uni- 
verse. 

The most direct way to constrain the star formation 
and chemical evolution history over cosmic times is to trace 
back galaxy properties (star formation rate, metallicity, stel- 
lar mass) through observations at different redshifts. Sev- 
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eral studies on the evolution of the rest-frame UV emission 
density of galaxies, converted into star formation or metal 
ejection rates, have converged into a picture in which the 
maximum of galaxy star formation activity occurs over the 
redshift range 1 < z < 2 and declines sharply from z = 1 
towards the prese nt l|Lillv et al.lll996l : IConnollv et al.lll997l : 
ICowie et al.lll99"9l ). While several recent studies have built 
a consistent picture of the decline in cosmic star formation 
rate from z ~1 to the present (e.g. iHopkins fc Beacom|[2006l . 
and references therein), more uncertain is the behaviour at 
redshift higher than 2, because of the poor understanding 
of the effect of dust on the SFR derived from UV Spec- 
tral Energy Distributions (SEP) of h i gh-redshift galaxie s 
jMadau et all 1 19961 : ISteidel et ai]|l999l : 1lvison et al.ll2002l ). 



2 A. Gallazzi, J. Brinchmann, S. Chariot, S.D.M. White 



A broad peak of high star formation over the redshift range 
z ^ 1 — 2 and then a rapid decHne toward the present are fea- 
tures predicte d (or repr oduced) also by both ch emical cvo- 
lution models (|Pei fc Fall 1995: Ed munds fc PhiUip ps 1997; 
|Pei et aP 1 19991) and by semian alytic models of galaxy for- 
mation (e.g. iBaugh et al.lll998l ). 

A complementary approach is to study the chemical 
and star formation history over cosmic times through the 
so-called 'fossil cosmology', i.e. determining the past history 
of the Universe from its present contents. This approach 
has benefited from large spectroscopic surveys in the local 
Unive rse, such as the 2 dF Galaxy Redshift Survey (2dF- 
GRS. IColless et al-llioOll) and the Sloan Digital Sky Survey 
(SPSS. lYork et al.ri2000l ). which provide detailed spectral 
information for hundreds of thousa nds of galaxies. Based 
on suc h surveys, iBaldrv et al.l l|2002h and iGlazebrook et al.] 
l|2003l) have constrained the cosmic star formation history 
(SFH) from the 'cosmic optical spectrum', which represents 
the average emission from all the objects in a representa- 
tive volume of the Universe and has the advantage of being 
fitted by simpler models of star form ation historie s than 
thos e needed for individ ual objects. iHeavens et al.l (|2004l ) 
and I Jimenez et al] (|2005|) have applied a d ata compression 
algorithm (MOPED. iHeavens et al] l2000t ) to extract the 
SFH of ~ 100, 000 SDSS DRl galaxies from their optical 
spectra. This work has been recently e xtended to th e SDSS 
DR3 (three times larger sample) by IPanter et al.l (|200d ). 
with improvements both in the data and in the modelling 
sides (see also Tojeiro et al. 2007, Ocvirk et al. 2006 and Cid 
Fernandes et al. 2007 for similar methods to recover stellar 
content and star formation histories from galaxy spectra). 
This allowed them to derive the cosmic SFH from the 'fossil 
record' and study it as a function of the present-day stellar 
mass of galaxies. 

The contribution to the global SFR by galaxies of dif- 
ferent mass is being studied n ot only in th e local Uni- 
verse but also at higher redshift. I Juneau et al.l Hooi) stud- 
ied the dependence of the cosmic SFH directly on the 
stellar mass at the epoch of observation over the red- 
shift range 0.8 < z < 2, based on a near-infrared se- 
lected sample from th e Gemini D eep Deep Survey (G DDS, 
[Abraham et al.ll2004 ). Similarly, iBundv et al] l|2005l ) have 
quantified the decrease with redshift of the mass limit above 
which star formation appears to be quenched, based on a 
sample of more than 8000 galaxies in the redshift range 
0.4 < z < 1.4 drawn fr om the DEEP2 Galaxy Redshift 
Survey (iDavis et al.ll2003l). These resu l ts con firm those pre- 
viously found by iBrinchmann fc EllisI (|2000l ) for < z < 1. 
While high- and intermediate-mass galaxies have transi- 
tioned to a quiescent phase of star formation by z ~ 1, less 
massive systems dominate the star formation rate density 
till the present epoch. It appears, though, that the global 
star formation rate has been declining since z '^l for all 
galaxies populations, at a rate which is in dependent of stel- 
lar mass, as shown bv lZheng et"al] (I2007D estimating SFRs 
of ~15000 COMBO-17 galaxies from UV and IR luminosi- 
ties and accounting for individually IR-undetectable galaxies 
IZheng ct al. 2006). 

An important consistency check for all these studies 
comes from the comparison of the density of stellar mass 
and of metals at different epochs expected from the cosmic 
star formation and chemical enrichment histories (i.e. the in- 



tegral of these histories) with those directly measured. The 
evolution of the global st ellar mass density out to z = 3 has 
been first determined by [Dickinson et aL I l|2003l ). In concor- 
dance with estimates of the cosmic SFH, their study sug- 
gests that the redshift range 1 < z < 2.5 is a critical epoch 
when galaxies are growing rapidly attaining their final stel- 
lar mass. 

Much effort has been put also in measuring the 
chemical composition of galaxies at different epochs, 
through optical nebular emissi o n lin es studies at z < 



(e.g. iK obulnicky fc Zaritsky 'l99^; 'Lilly et al.l l2003l : 



iKobulnickv fc Kcwlcy 2004; Ellison ct al. 200^), through 
Lyman-break and UV-selected star- f orming galax i es up 



to z 



Pettini et ahl l200ll: [Steidel et al.l |2004| : 



IShaplev et ail |2004| : lErb et all |2006| ). through quasar 



absorption- line systems, in particular Damped Ly-a Ab- 
sorbers (DLA) at any redshift (e.g. [Pettini et al.l 19941 : 
Lanzetta et"allll995l : iPettini et a"i]|l997l : iPeroux et alll2003l . 



20051 . bOOet TAll these studies have highlighted a shortfall of 



metals in observed galaxy populations with respect to expec- 
tations from the cosmic SFH, known as the 'mis sing metals' 
problem at redshift around 2.5 ( Pettini 2006: Bouche et al.l 
[2005). A large fraction could be hosted in a recentl y dis- 
covered population of sub-DLAs jPeroux et al.l [20051 ) . Not 
more than 30 — 40 percent seems to be in intergalactic 
medium (IGM), and prob ably < 35 per cent of metals is still 
'missing' from the census l|Bouche et al.i. 2 007„) . These metals 
are likely locked in the hot gas phase ( Ferrara et al.l lioosi : 
iDave fc Oppenheimeij |2007| ) . The present-day distribution 
of metals is still highly uncertain, because little is known 
about the chemical composition of the possibly dominant 
baryonic component, the warm hot intergalactic medium 
(WHIM). However the fraction of metals contained in galax- 
ies, in particular those locked into the stellar component, 
has increased from 2 ~ 2 to the present, and is probably 
comparable to the fraction of metals outside galaxies (e.g. 
iDunne et al]|2003l : ICalura fc Matteucci|[20o3 ). 

In this work we focus on the baryonic and metal content 
of the stellar component of galaxies in the local Universe. 
What is the total amount of metals and baryons locked up 
into stars by the present epoch? What is the resulting av- 
erage stellar metallicity of the present-day Universe? To ad- 
dress these questions we join together information about the 
stellar mass and chemical properties of present-day galaxies, 
supported by the large statistics provided by the SDSS. The 
sample we analyse span large ranges in physical, spectral 
and morphological properties, and constitute in this sense a 
representative sample of the local Universe. This allows us 
also to study how metals and baryons are distributed among 
galaxies with different properties. In particular, we want to 
quantify the fraction of metals, in comparison to the frac- 
tion of baryons, locked up in galaxies as a function of their 
stellar mass, morphology and age. 

We exploit new estimates of physical parameters, such 
as stella r metallicity and ste llar mass, that we previously 
derived (jGallazzi et al.ll2005l , hereafter paper I) for a large 
sample of nearly 2 x 10^ galaxies drawn from the Sloan Dig- 
ital Sky Survey Data Release Two (SDSS DR2). We include 
all galaxies types, from quiescent early-type to actively star- 
forming galaxies. In our previous works we focused only on 
galaxies with high signal-to-noise ratio (S/N) spectra, be- 
cause of the poor constraints that can be obtained on stellar 
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metallicity from low-S/N spectra. We circumvent here this 
problem by stacking individual spectra of low-S/N galaxies 
with similar properties in order to obtain high-S/N (aver- 
age) spectra. 

The sample analysed is described in Section [2.11 along 
with the stacking technique adopted in order to include 
galaxies with low-S/N spectra (Section 12. 2p . We illustrate 
the physical parameters estimates extracted from individ- 
ual galaxy spectra and from coadded spectra in Section [2.31 
In Section 13.11 we derive the mass density of baryons and 
of metals locked up in stars, expressed also in terms of the 
average stellar metallicity of the local Universe. We then 
discuss several sources of systematic uncertainties in Sec- 
tion |3]2] and compare with other estimates in the literature 
in Section l3^ S ect ion |4?T] provides an inventory of the stellar 
metallicity and stellar mass today, focusing on the charac- 
teristic age of the stellar mass and metallicity distributions 
in Section 14.21 We compare the observed distributions of 
stellar mass and stellar metallicity with those predicted by 
the Millennium Simulation in Section 14.31 We finally sum- 
marise and conclude in Section [5] Throughout the paper 
we adopt a flat cosmology with = 0.3, flA ~ 0.7 and 
Ho = 70 hro km s~^ Mpc~^. The mo dels used for this work 
are computed for a IChabrierl (|2003l ) initial mass function 
(IMF) and are based on a metallicity scale where the solar 
metallicity is Z© = 0.02. AH the magnitudes used in this 
work are SDSS model magnitudes, unless otherwise speci- 
fied. 



2 THE APPROACH 

In this section we give a brief overview of the sample anal- 
ysed and the method applied to derive estimates of physi- 
cal parameters, such as stellar metallicity, (light- and mass- 
weighted) age and stellar mass (Section 12.11 the reader is 
referred to paper I for a more thorough description of the 
method). The method requires spectra with high S/N. To 
derive a fair estimate of the total budget of mass and metals 
in stars today we need however to include all objects. We 
include low-S/N galaxies by adopting a stacking technique, 
which is described in Section 12.21 We compare the physi- 
cal parameters of the coadded spectra to those of individual 
galaxies in Section [2.31 

2.1 The sample 

To derive an estimate of the total amount of metals and 
baryons locked up in stars today and to study their distri- 
bution as a function of various galaxy properties we exploit 
a large sample of galaxies, for which stellar metallicities, 
as well as other physical parameters, have been estimated. 
The sample analysed here is dra wn from the main spe ctro- 
scopic sample of the SDSS DR2 (|Abazaiian et aLlliooi ) and 
is based on 164,746 unique spectra of galaxies with Petrosian 
r-band magnitudes in the range 14.5 ^ r ^ 17.77 (after cor- 
re ction for Galactic ex tinction using the extinction maps 
of ISchlegel et al.lll998l) . and with redshiftQ between 0.005 

^ As explained in paper I, we choose to limit the analysis in this 
redshift range, in order to avoid redshifts for which deviations 



and 0.22. The sample includes all galaxy types, from star- 
forming late-type to quiescent early-type galaxies. We note 
that the sample analysed is defined on the DR2 coverage, 
but we use the photometric reduction of the DR4 release. 
This is motivated by the fact that we found a systematic 
difference of ~0.16 mag in the z-band model magnitudes 
for a subset of galaxies from one release to the other, which 
can affect the overall normalization of the stellar mass den- 
sity. The spectroscopic measurements and fibre colours were 
instead consistent within the errors between releases0 

Bayesian-likelihood estimates of the stellar metallicities, 
r-band light-weighted ages and stellar masses of the galax- 
ies in the sample have been obtained in our previous work, 
by comparing the s pectr um of each galaxy to a library of 
iBruzual fc CharlotI (|2003l . hereafter BC03) models, covering 
the full range of physically plausible star formation histo- 
ries. The comparison is based on five spectral absorption 
features, namely D4000, H/3 and H5a+H7a as age-sensitive 
indices, and [Mg2Fe] and [MgFe]' as metal-sensitive indices, 
all of which have at most a weak dependence on element 
abundance ratios. After constructing the probability den- 
sity function of age, metallicity and stellar mass for every 
galaxy, the median of each likelihood distribution represents 
our estimate of the corresponding parameter, while half of 
the 16 — 84 percent interpercentile range gives the associated 
±1(7 (Gaussian-equivalent) uncertainty. In this work we add 
information about the mass- weighted age of galaxies. In Sec- 
tion 14.31 we shall use this quantity also in co mparison with 
predi ctions from the Millennium Simulation (|Springel et al.l 
l2005h . We have derived mass- weighted ages in the same way 
as the other physical parameters as described in paper I. The 
mass-weighted age of each model in the library has been es- 
timated by weighting each generation of stars by their mass, 
taking into account the fraction of mass returned to the in- 
terstell ar medium (ISM) by long- lived stars. 

In lGallazzi et al.l l|2005l . |2006| . hereafter paper II) we fo- 
cused only on galaxy spectra with median S/N per pixel 
of at least 20. As explained there, this is the minimum S/N 
required in order to obtain reliable estimates of stellar metal- 
licity. The quality of the spectrum infiuences directly the un- 
certainties in the derived physical parameters, stellar metal- 
licity being the most affected one: the average uncertainty 
on stellar metallicity decreases from 0.21 dex to 0.12 dex 
when high-S/N galaxies only are considered. The cut in S/N 
excludes roughly 75 percent of the galaxies and biases the 
sample towards high-surface brightness, high-concentration, 
low-redshift galaxies. Only 10 percent of the galaxies with 
concentration paramete 10 C ^ 2.4 satisfies the S/N require- 
ment. Excluding galaxies with S/N < 20 we would therefore 
preferentially miss diffuse systems with potentially subsolar 
metallicity. In order to derive a fair estimate of the total 
metal budget in the local Universe we need to include all 



from the Hubble flow can be substantial and to include galaxies 
in the stellar mass range 10® — 10^^ Mq with a signal-to-noise per 
pixel of at least 20. 

^ The stellar metallicity, light-weighted age and stel- 
lar mass estimates for the whole DR4 are available at 
http: //www. mpa-garching.mpg.de/SDSS/DR4/ 

defined as the ratio between the radii including 90 and 50 per- 
cent of the r-band Petrosian flux. 



4 A. Gallazzi, J. Brinchmann, S. Chariot, S.D.M. White 



galaxies down to the magnitude limit of the survey, there- 
fore low-S/N galaxies need to be considered as well. 

2.2 The stacking technique 

In order to include low-S/N galaxies, in addition to the sub- 
sample with S/N 20, we create composite high-S/N spec- 
tra by coadding the spectra of low-S/N galaxies with similar 
properties. First of all, we require galaxies to have similar 
velocity dispersion. The broadening due to stellar velocity 
dispersion affects the measured spectral absorption indices. 
When deriving physical parameters estimates we do not cor- 
rect for this, instead each spectrum is compared only to 
those models in the library with velocity dispersion similar 
to the observed one. It is therefore important that the galax- 
ies that contribute to each coadded spectrum span a range 
in velocity dispersion comparable to the observational er- 
ror. Moreover, metallicity, age and stellar mass all show cor- 
relations with velocity dispersion, absolute magnitude and 
D4000 (see e.g. figs 7,8 of paper I and figs 6,10 of paper II 
for early-type galaxies) . We thus choose to coadd the spectra 
of low-S/N galaxies with similar velocity dispersion, r-band 
absolute magnitude and 4000A-break. By binning into these 
quantities we are confident that the scatter in the physical 
parameters of the galaxies contributing to each stacked spec- 
trum is small. 

We first divide galaxies into bins of velocity dispersion 
logov of width /S.\ogav = 0.05 and bins of r-band absolute 
magnitude Mr of width AMr = 0.5. In each of such bins, 
galaxies are then ordered with increasing D4000 strength, 
and their spectra are stacked until a minimum S/N of 40 is 
reached. Each spectrum is weighted by 1/Vmax, where Vmax 
is the maximum visibility volume given by the bright and 
faint magnitude limits of the sample (14.5 ^ r ^ 17.77), and 
by our requirement that the galaxy redshift be included be- 
tween 0.005 and 0.22. The true number density of galaxies in 
the Universe should be estimated by accounting for galaxies 
that are missed due to, e.g., fibre collisions and spectroscopic 
failures. To correct for this, we have compared the r-band 
luminosity function obtai ned with our Vmax es timates with 
the luminosity function of lBlanton et al.l ()2003l ) and derived 
a normalisation factor for our Vmax estimates. At the end 
we obtain 14,694 coadded spectra from 122,643 spectra of 
low-S/N galaxies in the redshift interval 0.005 < z ^ 0.22. 

Fig. [T^,b show the distribution in velocity dispersion 
and r-band absolute magnitude for the coadded spectra 
(solid line), compared to the distribution for the individual 
low-S/N galaxies (dot-dashed line). For each stacked spec- 
trum we estimate the absolute magnitude Mj in a band j 
as the weighted sum of the luminosities Li^j of the low-S/N 
galaxies contributing to the coadded spectrum, according to: 

M, = -2.5 log + const. (1) 

where Wi is the weight 1/Vmax of the individual galaxies. 
The distribution in these two quantities as obtained from 
the stacked spectra agrees very well with the original distri- 
bution for the low-S/N galaxies, as expected since galax- 
ies have been binned in velocity dispersion and absolute 
magnitude. It is interesting to look how well the distri- 
bution in other morphological and photometric properties, 
into which galaxies are not explicitely binned, is reproduced. 



Fig. [T};,d show the distribution in the concentration param- 
eter C = -R90/-R50, where Rgo and -R50 are the r-band Pet- 
rosian radii, and in rest-frame g — r colour. The colour of 
each stacked spectrum is estimated as the difference between 
magnitudes defined according to equation [1] The concentra- 
tion parameter assigned to each stacked spectrum is given 
by the 1 /Vmax weighted-average concentration parameter of 
the galaxies that contribute to the stacked spectruirlf|. The 
distributions for the stacked spectra and the low-S/N galax- 
ies agree reasonably well, due to the correlation between 
colour and velocity dispersion or magnitude, and the small 
scatter in concentration parameter at given \ogav, Mr and 
D4000 (the mean absolute deviation in each such bin is typ- 
ically 0.18). The dotted line in each panel of Fig.[l]shows for 
comparison the distribution for the high-S/N galaxies. This 
clearly shows that by excluding low-S/N galaxies we would 
miss a substantial fraction of small, low-concentration, blue 
galaxies, i.e. preferentially young, metal-poor, star forming 
galaxies. 

We note that the distribution in concentration param- 
eter obtained from the coadded spectra is clearly bimodal 
and narrower than the distribution of the original low-S/N 
sample. A bimodality in C is expected from the bivariate 
distribution in the plane described by C versus D4000. The 
choice of stacking spectra with respect to D4000and the def- 
inition of weighted-average concentration parameter for the 
coadded spectra give higher-S/N measures of the concentra- 
tion index for 'blue' sequence and 'red' sequence galaxies 
separately, thus enhancing the bimodality in C. 

From each stacked spectrum, we also measure D4000, 
the higher-order Balmer lines and the other spectral absorp- 
tion indices defined in the Lick system, in the same way as 
they are measured from the spectrum of individual galax- 
ies (see also section 2.2 of paper I). They represent the 
1/ Vmax-weighted average of the absorption indices of the 
galaxies that contribute to each coadded spectrum. More 
properly, the fluxes in the 'pseudo-continuum' and central 
bandpasses measured from the coadded spectrum are the 
1/ Vmax-weighted average of the fluxes measured from the in- 
dividual galaxy spectra. In Fig.[2]the distribution in the flve 
spectral absorption features used to constrain stellar metal- 
licity, age and stellar mass estimates as measured from the 
stacked spectra (solid line) is compared to the distribution 
for the original sample of 122,677 low-S/N galaxies (dot- 
dashed line). The distributions for the stacked spectra are 
in very good agreement with the distributions for the orig- 
inal low-S/N galaxies. This is particularly true for D4000, 
as expected, since the spectral coaddition is performed on 
galaxies with similar D4000. The comparison for the other 
indices shows that the increased signal-to-noise ratio in the 
stacked spectra removes the tails of outliers present in the 
distributions for the original low-S/N galaxies, but absent 
in the distributions for the 42,103 high-S/N galaxies. 



* Very similar results are obtained if we assign to each stacked 
spectrum a concentration parameter given by tlie ratio between 
the weighted-average Petrosian radii. 
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Figure 2. Distribution in the five spectral absorption features adopted to derive estimates of stellar metallicity, age and stellar mass for 
the sample of 14,694 high-S/N stacked spectra (solid line), corresponding to 122,643 low-S/N galaxies (dot-dashed line). The absorption 
indices are measured off the stacked spectra in the same way as they are measured off the spectrum of individual galaxies and they 
represent the 1/Vmax-weighted average index of all the galaxies that enter each stacked spectrum. Dotted lines represent the distribution 
for the 42,103 high-S/N galaxies. 
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Figure 1. Distribution in velocity dispersion (a), r-band absolute 
magnitude (b), concentration parameter (c), and (fc-corrected) 
g—r colour (d). The solid line in each panel shows the distribution 
for the sample of 14,694 stacked spectra. This can be compared 
to the distribution of the low-S/N galaxies, shown by the dot- 
dashed line. The dotted line in each panel represents instead the 
distribution of the high-S/N galaxies. 

2.3 Physical parameters estimates 

Estimates of stellar metallicity, (light- and mass-weighted) 
age and stellar mass are derived from the coadded spectra 
in the same way as they are derived from individual galaxy 
spectra, as summarised in Section [2. II and more extensively 
described in paper I. The physical parameters are derived 



by fitting the galaxy spectra as observed and so they re- 
fer to the galaxies at the time of observation. This concerns 
in particular the stellar age. In Section 14.11 and 14.21 we will 
study the distribution of metals as a function of stellar age. 
In order to define a characteristic age and interpret it as a 
characteristic redshift of metal production, we correct the 
measured mass- and light- weighted ages by adding the look- 
back time to the redshift at which the galaxy is observed. 
The age obtained in this way represents the effective (mass- 
or light-weighted) epoch when stars formed. For the stacked 
spectra we assume the average redshift of all the galaxies 
that contribute to each spectrum. The spread in redshift of 
the galaxies contributing to each coadded spectrum is on 
average 30 percent for redshift up to 0.1 and 20 percent for 
z > 0.1. 

The distribution in the derived parameters for the whole 
sample of 164,746 galaxies is shown in the left-hand pan- 
els of Fig. [3] (thick solid line). The mass- weighted age of 
the sample peaks at ~10 Gyr and then extends towards 
ages as young as 2.5 Gyr. The r-band light- weighted age 
shows a roughly bimodal distribution with a primary peak 
around 9 Gyr and a broader peak around 4 Gyr. The dis- 
tribution in stellar metallicity is also highly skewed with 
a primary peak around 1.4 x Zq and an extended tail to- 
wards lower metallicities. The effect is much weaker for stel- 
lar mass, which has a distribution roughly symmetric around 
a mean log(M«/MQ) = 10.81 with a scatter of 0.46 dex. We 
note that, as expected, the derived mass-weighted age of 
the galaxies is older than their light-weighted age. This dif- 
ference is larger for younger galaxies, going from 0.7 Gyr 
for galaxies with tr=6.3 Gyr up to 4 Gyr for galaxies with 
tr-=2.5 Gyr. This reflects (on average) the more extended 
star formation histories of younger (less massive) galaxies. 

The right-hand panels show the distribution in the un- 
certainties on mass- weighted age, r-band light-weighted age, 
stellar metallicity and stellar mass, given by half of the 
16 — 84 percent percentile range of the corresponding like- 
lihood distribution. The dotted line shows the distribution 
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Figure 3. Distribution in mass-weiglited age, r-band light- 
weighted age, stellar metallicity and stellar mass (from top to 
bottom, left panels) for the final sample obtained by combining 
high-S/N galaxies and the low-S/N galaxies included in the coad- 
ded spectra (thick solid line). The dotted lines show the contribu- 
tion by high-S/N galaxies only, while the dot-dashed lines repre- 
sent the distribution in the parameters for the low-S/N galaxies, 
as derived from the stacked spectra. The right-hand panels show 
the distribution in the corresponding uncertainties, given by half 
of the 16 — 84 percent percentile range of the likelihood distribu- 
tion. The grey-shaded histograms in the right-hand panels give 
for comparison the distribution in the 68 percent confidence in- 
terval of the physical parameters of low-S/N galaxies as derived 
from their individual spectrum. 

for the high-S/N galaxies, while the dot-dashed line shows 
the distribution for the uncertainties on the parameters of 
low-S /N galaxies as derived from the coadded spectra. This 
can be compared to the 68 percent confidence range on the 
parameters of low-S/N galaxies as derived from their indi- 
vidual spectra (grey-shaded histogram). This makes clear 
the importance of a good spectral S /N in the determination 
of the physical parameters (in particular stellar metallicity, 
see also paper I) and the advantage of the stacking tech- 
nique: it allows us to retrieve the physical parameters of 
galaxies with low-S /N spectra with a much better accuracy 
than what we could do from their individual spectra. 

The physical parameters derived from the stacked 
spectra can be interpreted as the (l/Vmax)-weighted aver- 
age stellar metallicity, age and stellar mass of the galaxies 
that contribute to each coadded spectrum. To test how 
well we can recover the physical parameters of individual 
galaxies with our stacking technique, we have generated a 
control sample of stacked spectra by coadding the spectra 
of individual high-S /N galaxies, for which reliable estimates 
of metallicity, age and mass can be derived, in the same 



way as described in Section 12.21 for the low-S /N galaxies. 
Before coaddition we added gaussian noise to the individual 
high-S/N spectra to mimic the situation we have when 
coadding low-S/N spectra. We then compare the physical 
parameters estimated from the coadded spectra with the 
(l/Vmax)-weighted average parameters of the galaxies 
that contribute to each coadded spectrum. This is shown 
in Fig. |4] for stellar metallicity (panel a), stellar mass 
(panel b), light- weighted age (panel c) and mass- weighted 
age (panel d). The histogram of the difference between 
the derived ('stack') and expected ('wavg') parameters 
is compared to gaussian distributions of width given by 
the average uncertainty on the derived parameter (dashed 
line) and by the average scatter in the physical parameter 
of the galaxies that contribute to each stacked spectrum 
(dot-dashed line). For all the parameters we can recover 
the expected value within the corresponding typical error. 
We note however that there is a small but systematic offset 
(dotted vertical line) on average of about —0.02 dex in 
stellar mass, —0.03 dex in mass-weighted age and —0.04 dex 
in light-weighted age (while the average offset in stellar 
metallicity is negligible). This offset likely originates from 
the fact that each coadded spectrum tend to be dominated 
by the youngest, brightest stellar populations of the galaxies 
that compose it, and therefore the derived light-weighted 
ages and mass-to-light ratios tend to be biased low. We take 
this into account as a systematic uncertainty, as discussed 
in Section [3.21 below. 



3 THE MASS DENSITY OF BARYONS AND 
METALS IN THE LOCAL UNIVERSE 

In this section we derive an estimate of the total amount of 
baryons and metals locked up in stars and the average stellar 
metallicity in the local Universe, by combining the contribu- 
tion of individual high-S/N galaxies and the low-S/N galax- 
ies included in the coadded high-S/N spectra (Section 13. If) . 
We also discuss and quantify systematic uncertainties in Sec- 
tion [3]2] and compare with various observational estimates 
and model predictions in the literature in Section 13.31 

3.1 The total stellar metallicity in the local 
Universe 

We compute the mass density of metals in stars, pz, and 
the stellar mass density, p*, at a mean redshift oi z = 0.1, 
as follows: 

pz = Y. ^-'i + E wt) (2) 

i % 

p,^J2 + E ^^') (3) 

i i 

where Z« and M« are the median-likelihood estimates of the 
stellar metallicity and stellar mass of the high-S/N galaxies, 
and w are the weights 1/Vmax. The symbols Zf and Mf 

^ At low S/N there is likely to be non-gaussian noise sources 
such as sky subtraction problems, but our approach here should 
capture most of the trends. 
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Figure 4. Distribution of the difference between the stellax metal- 
licity (a), stellar mass (b), light- weighted age (c) and mass- 
weighted age (d) estimated from the stacked spectra and the 
(l/Vmax)-weighted average parameter (indicated with 'wavg') 
of the galaxies that contribute to each stacked spectrum. Each 
histogram is compared to gaussian distributions of width given 
by the average uncertainty on the corresponding parameter ((cr), 
dashed line) and by the average rms scatter in the physical pa- 
rameters of the galaxies that enter each stacked spectrum (rms, 
dot-dashed line). For this test we used a control sample of stacked 
spectra obtained by coadding the spectra of individual high-S/N 
galaxies, for which reliable estimates of metallicity, mass and age 
can be obtained, after degrading their quality adding gaussian 
noise. 



refer to the stellar metallicity and mass estimated from the 
coadded spectra. These have to be weighted by W"*, which 
is the sum of all the weights 1/Vmax of the low-S/N galaxies 
contributing to each stacked spectrum. 

Combining the contribution of individual high-S/N 
galaxies and coadded spectra, we derive: 

pz = 7.099 ± O.OWttllt X 10*^ h7o M© Mpc"^ (4) 

p. = 3.413 ± 0.005to.554 X 10** /iro Mq Mpc"^ (5) 

We note that the derived stellar mass density corresponds 
to 0.0025 times the critical density (see Table ??). We thus 
find a stellar baryon fraction of only 7 percent, assuming a 
value of 0.004 for the cosmic baryon fraction. 

The systematic uncertainties on these quantities (and 
on the average stellar metallicity {Z,/Zq) derived below) are 
discussed in detail in Section|3]2]and summarized in Table|3l 
We calculate 'corrected' mass densities using stellar mass 
and stellar metallicity estimates corrected for systematics as 
described in Section |3]2l The uncertainties on p* and pz are 
then expressed as the difference between the 'corrected' and 
the default values. The sources of systematic uncertainties 
that we consider include the possible bias introduced by the 
adopted stacking technique, the non-modeled dependence 
of index strengths on element abundance ratios, the prior 
distribution of the models in the Monte Carlo library, the 



aperture effects and the magnitude used to normalize the 
stellar mass-to-light ratio. 

The systematic uncertainties quoted in equations |3] 
and [S] are the largest among those summarized in Table O 
These are the aperture bias, which can lead to an overes- 
timate of ~16 percent and of ~27 percent in p« and pz 
respectively, and the stacking technique, which may intro- 
duce an underestimate in pz of about 30 percent and less 
than 20 percent in p, . The use of Petrosian magnitudes for 
determining stellar masses would provide estimates of both 
p* and pz lower by less than 13 percent of those derived 
with model magnitudes. A similar effect would be caused 
by the choice of a burst-enhanced prior. For completeness in 
equations |4] and [5] we also indicate the statistical error (esti- 
mated by standard error propagation from the uncertainties 
on metallicity and mass). This is very small and represents 
only the 0.2 percent of the value of pz and 0.1 percent of 
the value of p« that we derive. 

We can now combine Equation [5] and [3] to estimate the 
(mass-weighted) average metallicity in stars in the nearby 
Universe, obtaining: 



{Z^ M.,. w,) + {zf M:^i Wj') 
1.04±0.002lo;if ^0 



(6) 



The normalization of the total metallicity in stars is affected 
mainly by the choice of prior and by aperture effects, which 
can lead to an overestimate and underestimate, respectively, 
of about 14 percent. The uncertainties introduced by the 
stacking technique amount to roughly 12 percent (see Ta- 
ble |3J. The scaled-solar abundance ratio of the models and 
the choice of normalization magnitude can introduce sys- 
tematics lower than 10 percent. We note that, when consid- 
ered separately, high-S/N galaxies give a total metallicity of 
1.149^0 while low-S/N galaxies give, as expected, a lower 
metallicity of 0.996^0 

We derive thus that the average metallicity of stars in 
the present-day Universe is the typical metallicity of L* 
galaxies, i.e. consistent with the solar value. We mention 
in passing that, while the solar metallicity scale has been 
substantially adjusted down wards with respect to the pre- 
vious ly recommended value l|Asplund et al.l l2005l : lAsplundl 
l2005h . this has no impact on our analysis. This is because 
the BC03 models are tied to the iron abundance [Fe/H], 
whose solar value is unchanged with the new calibration. In 
view of this we adopt Zq = 0.02 in this paper for consistency 
with our previous work. 

Our result, based on a large homogeneous sample of 
galaxies spanning large ranges in physical properties, and 
on recent population synthesis models accounting for the 
full range of possible SFHs, puts on a more robust basis 
(and with accurate estimates of the systematic uncertain- 



^ Low-S/N galaxies are preferentially galaxies with low surface 
brightness and low concentration parameters. As shown in pa- 
per I, while (the bulges of) massive disc-dominated galaxies have 
metallicities comparable to those of bulge-dominated galaxies, 
low-mass C ^ 2.4 galaxies are on average more metal-poor than 
their C 2.8 counterparts. We thus expect low-S/N galaxies to 
have a mass-weighted mean metallicity slightly lower than high- 
S/N (mostly bulge-dominated) galaxies. 
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ties) an early calculation by lEdmunds fc Phillippd |l993), 
based on closed-box chemical models for irregulars, ellipti- 
cals and spirals. Assuming mean metal abundances and com- 
bining contribution by diflterent galaxy types according to 
their different luminosity functions and average stellar mass- 
lo- light ratios, they derive a mass- weighted mean oxygen 
abundance (including stars and gas) of 12-|-log(0/H) = 8.8, 
i.e. consistent with s olaj^ . This is also in agreement with 
ICalura fc Matteuccil (|2004 '). who, summing the contribution 
of ellipticals, spirals and irregular galaxies from a chemo- 
photometric model, predict a mass-weighted mean metallic- 
ity in stars of 0.019, i.e. roughly 0.95 Zq (from their tables 
1, 2, 3, 9). 

We stress that the metal mass density that we estimate 
here accounts only for the fraction of metals in the stellar 
populations of present-day galaxies, and not those in the in- 
terstellar medium (nor in other gaseous form outside galax- 
ies). In the local Universe, the stellar components in galaxies 
host a significant fraction (30 — 50 percent) of the cosmic 
metal content, almost comparable to the fraction of metals 
that reside outside galaxies (according to Calura&Matteucci 
2004, and as recently confirmed, within uncertainties, by 
Dave&Oppenheimer 2007 using cosmological hydrodynam- 
ical simulations which incorporates enriched galactic out- 
flows in a hierarchical structure formation scenario). This 
is not the case at higher redshift (above 2 ~ 2), when the 
diffuse inter-galactic gas is the largest reservoir of metals. 
These considera tions find support in the observational cen- 
sus compiled bv lDunne et al.l l|2003l ). 

Following naiv ely the discussion of 

lEdmunds fc PhilhppsI (|l997l ) the solar average stellar 
metallicity in the local Universe would be consistent 
with the idea that most of the baryons available for star 
formation have been locked up into stars (if the effective 
yield has approximately solar value) Considering only 
galactic components, while at high redshift the majority 
of baryons are in the ISM, at low redshif t the largest 
reservoir of baryons are the s tars (see e.g. IPunne et al.l 
I2OO3I : ICalura fc Matteuccilliooi ) . This is consistent with the 
decline in the ensemble star formation rate since 2: ~ 1 till 
the present. However, the majority of all baryons appears 
to be at all epochs outside g alaxies, in the diffuse IGM (e.g. 
iDave fc Qppenheimerir2007h . 



3.2 Sources of uncertainties 

We discuss here the possible sources of systematic un- 
certainties that we identify in our procedure for deriving 
physical parameters estimates from individual or coadded 
galaxy spectra (Tables [1] and [2}, and quantify how they 
propagate into the estimates of stellar mass density, mass 
density of metals in stars, and average metallicity (Table|3]). 
The results quoted in Equations 14] O |6] are not corrected 
for systematics. For each source of uncertainty, we esti- 
mate here, as best as we can, corrections on individual 



This value would correspond to 1.28 times the solar value, 
adopting a solar oxygen abundance of 12 + log(O/H)0 = 8.69 
according to lAUende Prieto et al.l l|200lh . 

* In a closed-box system the average stellar metallicity would be 
approaching the yield as the gas is progressively processed into 
stars. 



stellar masses and metallicities (and ages), and calculate 
'corrected' p, , pz and (2',). The systematic uncertainties 
on the these quantities are then expressed as the difference 
between their value after correction and the default one (as 
given in Section |3TT]) . Fig. |9]in Section \AJ2\ provides a visual 
representation of the extent of each systematic uncertainty 
as a function of galaxy stellar mass. 

Stacking technique. We already mentioned that the 
stacking technique we adopt may provide physical param- 
eters estimates which are (slightly) systematically lower 
(on average) than the expected average parameters of 
low-S/N galaxies. This affects in particular stellar mass 
and mean ages. To quantify the extent of these effects we 
measured the range between the upper and lower quartiles 
of the distributions shown in Fig. |4l i.e. the difference 
between measured and expected parameters (first row of 
Table [TJ. We do this also for stellar metallicity although the 
average offset is negligible. As mentioned before, while the 
differences in stellar metallicity are symmetric around zero 
and well within the typical uncertainty on metallicity, the 
masses and ages estimated from the coadded spectra are in 
general lower than the expected parameters by an amount 
corresponding to ~30 percent of the typical uncertainties 
on these parameters. We correct stellar metallicity and 
mass estimates according to these offsets and calculate 
'corrected' p*, pz and (Z,). The systematic uncertainties 
on these quantities are then expressed as the (interquartile 
range of the) difference between the 'corrected' and default 
values (Table [31. 

a/Fe abundance ratio. The Monte Carlo library we use 
is based on the BC03 population synthesis code, which does 
not include the effect of variations in element abundance 
ratios with respect to solar. In deriving galaxy physical pa- 
rameters we take care to constrain only those absorption in- 
dices which have a weak dependence on a/Fe abundance ra- 
tio. Nevertheless, recent studies have shown that the higher- 
order Balmer lines (also used in our analysis) are sensitive 
to a-enhancement at metallicities around solar and above 
and may lead to an underestimate of stellar ages if models 
with scaled-solar abundance ratios are used (Frii omas et al.l 
I2OO4I : iKorn et aLllioOsI : iProchaska erallbOOTl ). We estimate 
the possible (mass-dependent) systematic uncertainty on the 
physical parameters we derive with scaled-solar abundance 
ratios models as follows. 

First of all, for a stellar population with abundance ratio 
[a/Fe] = 0.3, we determine the offset between the expected 
(true) parameters and those derived with a model library 
with [a/Fe] = 0, using the predictions from Irhomas et al.l 
(|2004h simple stellar populations (SSP) models. This is de- 
scribed in detail in section 2.4.2 of paper I. The second row 
of Table [T] gives the interquartile range of the offsets (for 
[a/Fe] = 0.3) in stellar metallicity, light- and mass-weighted 
age, and stellar mass-to-light ratio. 

The degree of a-enhancement in ellipticals is a function 
of the galaxy stellar mass. Therefore we derive mass- 
dependent corrections to the physical parameters on a 
galaxy- by- galaxy basis and quantify the effect on p», pz 
and (Zt). We cannot derive an estimate of the true a/Fe 
abundance ratio of individual galaxies with our models. We 
use instead an empirical estimate given by A(Mgb/(Fe)), 
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i.e. the difference between the observed Mgb/{Fe) index 
ratio and the index ratio of the [a/Fe] — best-fit model 
in our hbrary (see section 2 of paper II). We then use 
the relationship between A(Mgb/(Fe)) and stellar mass 
derived for bulge-dominated galaxies (C ^ 2.8, see table 4 
of paper II). This gives us the degree of a/Fe expressed in 
terms of the galaxy stellar mass. From this we estimate the 
correction to be applied to the galaxy physical parameters, 
scaling the offsets derived above (Table [TJ accordingly. We 
apply these corrections only to C ^ 2.8 galaxies. 

Prior. Another possible source of systematic error 
comes from the choice of prior according to which our 
model library populates the parameter space, in particular 
the mixture of bursty and continuous star formation histo- 
ries. As discussed in section 2.4.2 of paper I, we explored 
the effect on the derived physical parameters of changing 
the fraction of bursts in the last 2 Gyr from 10 percent 
(our standard prior) to 50 percent. The mean ages and the 
stellar mass derived with this burst- enhanced prior are on 
average lower than those derived with our standard prior, 
while the stellar metallicity is slightly higher. We quantify 
this as the interquartile range of the difference between 
our standard prior and the burst-enhanced one in stellar 
metallicity, light- and mass-weighted age, and stellar mass 
(third row of Table As above, we obtain the range of 
the corresponding uncertainty on p,, pz and (Z*) as the 
difference between the 'corrected' and default values. 

Aperture effects. Estimates of stellar metallicity and 
mean stellar age are affected by the aperture bias, due to 
the fact that the SDSS spectra sample only a limited inner 
region of the galaxy. The light collected by the fibre is on 
average 30 percent of the total flux, but this fraction de- 
pends on the stellar mass, morphology and redshift of the 
galaxy. Due to the presence of stellar populations gradients, 
the metallicities and ages derived from the SDSS spectra are 
not representative of the galaxy as a whole but only of the 
bulge or central regions. A correction for this would require 
an accurate knowledge of metallicity gradients as a function 
of galaxy type and mass, and is not feasible here. This is 
clearly a concern in this work, where we want to estimate 
the total amount of metals in stars today (see also section 
3.4 of paper I for a discussion on aperture effects). 

We attempt to quantify this bias for galaxies with simi- 
lar stellar mass and concentration index by looking at trends 
in their estimated physical parameters as a function of the 
fraction of light in the fibre (given by the ratio / between 
the fibre and the Petrosian fiux). Stellar mass estimates can 
also be affected by aperture bias, because they are based on 
the (potentially wrong) assumption that the stellar mass- 
to-light ratio outside the fibre is the same as the one in- 
side (derived from the galaxy spectrum). In the same stellar 
mass and concentration bins as above we checked and quan- 
tified also the trend in the z-band stellar mass-to-light ratio 
log{Mt/ Lz) with /Ifl The linear relations fitted are summa- 
rized in Table [5] for those stellar mass and concentration 

^ For reference, a flux ratio of about 50 percent occurs when the 
fibre radius coincides with the effective Petrosian radius R50 for 
a de Vaucouleur proflle or is 1.3 X R50 for an exponential pro- 
file. When the fibre covers roughly twice R50, the flux gathered 



bins where the trends are statistically significant. The rela- 
tively large variation in metallicity that we find for bulge- 
dominated galaxies are, at least qualitatively, consistent 
with the typical abundance gradients fou nd in the literature 
. 16 — 0.2dex/decade in radius, e.g. iHenrv fc WorthevI 
ll999l : lMehlert et al.|[2003h . We note that we are not able to 
identify significant trends in stellar metallicity for C < 2.4 
galaxies. It is possible that we are underestimating the aper- 
ture effects for this class of galaxies. However, while abun- 
dance gradients are observed in individual spiral galaxies, 
they are on average rather shallow and associate d instead to 
large r age gradients (|MacArthur et al] f2004: Bel l fc de Jond 
I2OO j ). consistent with our findings. 

For each galaxy, we calculate new physical parameters 
estimates by adding the variation in the corresponding 
parameter obtained by extrapolating the relations given 
in Table [2] to / = 1. We then recalculate the total p., 
pz and {Z,}. As above, the systematic uncertainty on 
these quantities is expressed as the difference between the 
'corrected' and default values. 

Total magnitude. The stellar mass estimates are influ- 
enced by the total magnitude we choose as normalization 
of the stellar mass-to-light ratio. We use the SDSS model 
magnitudJ^ as a good estimate of the total galaxy lumi- 
nosity, with respect t o the more often adopt ed Petrosian 
magnitude (e.g. by iKauffmann et ahl |2003| ) . measured 
within a circular aperture of radius two times the Petrosian 
radius. The Petrosian flux should recover almost all the 
flux for an exponential profile and ~ 80 percent of the 
flux for a de Vaucouleurs proflle. For completeness and 
for comparison with other works, we also derive stellar 
mass estimates using the z-band Petrosian luminosity to 
normalize the M,/Lz. The last row of Table [1] gives the 
interquartile range of the difference between 'model' and 
'Petrosian' stellar masses. 

Finally, we note that the adopted shape of the initial 
mass function (IMF) and the assumption that it is univer- 
sally applicable between galaxies with different mass and 
star formation history clearly have an impact on the derived 
stellar mass, and hence p« and pz- Exploring in details the 
effects of the different assumptions goes beyond the scope 
of the present work. However, we try here to quantify how 
much the mass-to-light ratio scale would vary by varying 
the parameters of the IMF. In particular we consider the ef- 
fect of changing the Chabrier IMF parameters wit hin their 
quoted uncertainties (see Table 1 of IChabrierll2003l ) . We cal- 
culated the evolution of the z-band mass-to-light ratio for 
SSPs of different metallicity (from 0.4 to 2.5 times solar) for 
two 'modified' Chabrier IMFs, one with the steeper high- 
mass slope of 1.6 and upper limits for the low-mass end 
parameters, and another one with the shallower high-mass 
slope of 1 and lower limits for the low-mass end parameters. 
In the first case the ratio of Mt/L^ between the modified 
and the standard Chabrier IMF varies from 0.73 at 1 Gyr to 

is about 80 percent and 60 percent for a de Vaucouleur and ex- 
ponential proflle respectively. 

This is the total magnitude of the best-fit model between a 
de Vaucouleurs and an exponential proflle to the r-band image of 
the galaxy. 



10 A. Gallazzi, J. Brinchmann, S. Chariot, S.D.M. White 



0.66 at 10 Gyr. In the second case this ratio varies from 1.55 
to 1.75 over the same time interval. In both cases there is 
negligible dependence on metallicity. We also note that the 
difference in M^/Lz between the various IMFs roughly cor- 
responds to the difference between SSPs of different metal- 
licities at a given IMF. Similar results are obtained when 
looking at fixed absorption index strength, such as D4000 
or [MgFe]' or H(5,4+H7a, rather than at fixed age. We note 
that the the uncertainty in the slope at the high-mass end 
is largely responsible of the significant variation in M*/Lz 
between the 'modified' and 'standard' IMFs, by varying the 
weight of stars just above the turnoff mass of IMq . If we let 
vary only the low-mass end parameters, fixing the high-mass 
slope at 1.3, the M^/Lz is ^ 1.03 times and ^ 0.97 times (for 
'lower' and 'upper' case respectively) the M^jL^ of the stan- 
dard Chabrier IMF. Since it will be useful in Section[3]3] we 
also men t ion he re the comparison with the single power-law 
ISalpeteij ()l955l ) IMF. For a solar metallicity simple stellar 
population the ratio between the z-band h'h/L predicted 
by the Salpeter IMF and the one predicted by the Chabrier 
IMF is on average 1.75 and varies from 1.72 at 1 Gyr to 1.76 
at 10 Gyr. 



3.3 Comparison with the literature 

We list in Tables |4] an(45] the present-day mass densities of 
baryons and metals in stars estimated in this work together 
with the values derived from several sources in the litera- 
ture, distinguishing observational estimates and prediction 
from models or integration of the cosmic star formation his- 
tory. All the values in the literature have been converted into 
our adopted cosmology and to a Chabrier IMF, when neces- 
sarj0. The comparison is also illustrated in Figure (5] Over- 
all we find good agreement within our uncertainties with 
other observational estimates and model predictions. 

As far as the stellar mass density is concerned, 
our determinati o n is in line with the values derived by 
iKochanek et al.l (|200lf ) from the _ft'-band lumin osity func- 
tion o f a sample of ~5000 2MASS galaxies, by I Cole et al] 
lj200ll ) from the NIR luminosity function of a 2MASS and 
2dFGRS co mbined sample of infrared-selected galaxies, and 
similarly by iBell et all (|2003l ) from integration of the stel- 
lar mass function derived from optical-to-NIR SED fitting 
on a matched sample of SDSS/2MASS galax ies. We find 
good agreement also with iPanter et all (|2006l ). who apply 
the MOPED algorithm to fit the entire optical spectrum of 
SDSS DR3 galaxies, based on BC03 models. This is encour- 
aging given the similarity in the approach and in sample 
definition. The offset of ~0.07 dex is accounted for by the 
different definition of 'total' magnitude adopted (i.e. z-band 
model versus r-band Petrosian magnitude). 

Analysis based on the Millennium Galaxy Catalogue 
also provides stellar mass density e stimat es in agreement 
with ours, as shown in [Driver et al.l (|2007al ) . Recently, they 



FoIIowine lBell et al.l ||2003|) we increased by 0.15 dex the mass 
densities derived with a 'diet'-Salpeter IM F (i.e. a Salpeter IMF 
modified to have a flat slope below 0.6 M© iBeU fc de Jonidl200lf ) 
to convert them into a Salpeter IMF. We then assumed Salpeter- 
based stellar mass densities to be 1.75 times higher than Chabrier- 
based densities (see last paragraph of Section l3.2l l. 



revised their determination, by incorporating effects of inter- 
nal dust extinction based on an empirical relation between 
internal attenuatio n and inclin ation in galaxy discs and their 
associated bulges (iDriver et al. 2007h). This correction in- 
creases by ~20 percent their previously derived stellar mass 
density. Our estimates of p* are still in agreement within the 
combined uncertai nties. Similar considera tions hold for the 
value obtained by iFukugita et al.l (|l998h combining infor- 
mation on the luminosity densities and stellar mass-to-light 
ratios of spheroids, discs and irregular galaxies. 

iRudnick et al] l|2003D . based on a sample of SDSS EDR 
luminous galaxies, obtain instead a stellar mass density 
roughly 50 percent lower than our value; this discrepancy 
is only slightly alleviated with their more r ecent estimate 
accou nting for correction to 'total' values (jRudnick et al.l 
I2OO6I ). Finally, our estimate of stellar mass density falls in 
the range constrained by Glazcbrook et al. (2003), by fitting 
the SDSS-based 'cosmic optical spectrum' with a power-law 
SFH up to z = 1. 

We find marginal agreement also with the v alue o f 
stellar mass density estimated by IShankar et al.l l|2006l ). 
Their mass-to-light ratio estimates are based on a kine- 
matic decomposition for late-type galaxies (following 
ISalucci fc Persidll999|l and on central velo city dispersion for 
early- type galaxies i Bernardi et al.|[2003b[ ). without any as- 
sumption on IMF. 

Comparing our estimate of p* to model predictions, 
we find reason ably good agreement w ith the range of val- 
ues found by iNagamine et al.l (|2006l ). who explore four 
different approaches: a two-component Fossil model, an 
Eulerian hydrodynam i c sim ulation, the analytic SFH of 
iHernauist fc Springell (|2003l ) based on their cosmological 
smoothed particle h ydrodynam i cs sim ulations, and the semi- 
ajialytic model of ICole et al] l|2000t ). An early result by 
Pei et al.i (I1999h . based on a set of coupled equations linking 
the evolution of the densities of stars, gas, heavy elements 
and dust based on data from quasar absorption-line surveys, 
optical imaging and redshift surveys available at that time, 
is also in agreement within la with our estimate. We find 
instead a discrepancy of about 40 per cent with the stellar 
mass density derived by the models of lCalura fc Matteuccil 
(l2004h. A valu e lower by about 30 percent is presented by 
iMonaco et all (2007) based on the MORGANA code for the 
formation and evolution of galaxies and AGN. Finally, we 
compare our measured stellar mass density with the result 
of integration of the cosmic star formation history. The in- 
tegration (accounting for the fraction of mass r eturned to 
the IS M by evolved star s) of the analytic fi t of ICole et al-l 
lj200l[ ) to the SFH data of lSteidel et"al] (|l999t ) clearly under- 
estimate the total stellar mass density at the present epoch. 
On the other hand, integration of the dust-corrected SFH 
provides too high a stellar mass density (up to 70 percent 
higher). This discrepancy might be r elated to the sp ecific 
dust-correction applied. We note that iBell et al.l(|2007l l. ac- 
counting for both unobscured (from UV) and dust-obscured 
(from IR) star formation, find agreement between the inte- 
gral of the SFH and the measured stellar mass density. In 
Section 14.31 we will compare in detail our results with the 
predictions from the Millennium Simulation. We anticipate 
here that we find good agreement between the SDSS-derived 
stellar mass density and the stellar mass density predicted 
by the models, which amounts to 3.203 x lO^M© Mpc"^. 
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Table 1. Summary of the uncertainties or parameter choices that can affect the stellar metallicity (column 2), r-band light- weighted age 
(column 3), mass- weighted age (column 4) and stellar mass (column 5) estimates. The different entries, indicated in the first column, 
are: the offset between the parameters derived from the stacked spectra and the expected 1/Vmax-weighted average parameters of the 
corresponding low-S/N spectra; sensitivity to non-solar element abundance ratios, not taken into account in BC03 models; the prior 
according to which the models of the Monte Carlo library are distributed in parameter space (in particular we consider a library in which 
50 percent of the models had a burst in the last 2 Gyr against our default library in which this fraction is only 10 percent); the total 
2-band galaxy magnitude used to normalize the stellar mass-to-light ratio of the models (in particular we consider the difference of using 
Petrosian magnitude instead of our default choice of using model magnitude) . The extent of these uncertainties is expressed as the upper 
and lower quartiles of the difference between the measured value and the expected (corrected) one. 





A log(Z. /Zq) 


A\og{tr/yr) 


Alog{t M/yr) 


Alog(M,/MQ) 


(1) 


(2) 


(3) 


(4) 


(5) 


stacking 


-0.070 : 0.036 


-0.108 : 0.000 


-0.099 : 0.011 


-0.091 : -0.011 


o/Fe 


-0.084 : 0.146 


-0.166 : -0.003 


-0.170 : 0.026 


-0.119 : -0.014 


prior 


-0.057 : -0.010 


0.059 : 0.103 


0.051 : 0.091 


0.033 : 0.067 


total magnitude 








0.019 : 0.071 



Table 2. We quantify here the potential bias due to the fixed fibre aperture. We express this as the variation in each physical parameter 
as a function of the fraction of light missed by the fibre (A/) at given stellar mass and concentration parameter (column 1). In the table 
we give the slopes of those relations which are statistically significant. (See text for more details.) 



Range in C and log M* Alog(Z*/Z0) Alog(tr/j/r) A log(iA//s/r) AlogilVh/L^) 



(1) 








(2) 


(3) 






(4) 






(5) 


C ^ 2.8 


& 


10 < log(M,/M0) ^ 10.3 














-0 


32 ±0.13 


A/ 


C ^ 2.8 


& 


10.3 < log(M*/MQ) ^ 10.5 


-0.08 ±0.08 


A/ 










-0 


21 ±0.02 


A/ 


C > 2.8 


& 


10.5 < log(M*/MQ) ^ 10.8 


-0.14 ±0.06 


A/ 










-0 


13 ±0.02 


A/ 


C > 2.8 


& 


10.8 < log(M,/MQ) ^ 11 


-0.14 ±0.02 


A/ 










-0 


20 ±0.03 


A/ 


C ^ 2.8 


& 


log(M*/M0) > 11 


-0.19 ±0.03 


A/ 


-0.14 ±0.03 A/ 


-0.11 ±0 


02 


A/ 


-0 


16 ±0.02 


A/ 


2.4 < C 


< 


2.8 & 10 < log(M./M0) <: 10.3 






0.17 ±0.12 A/ 














2.4 < C 


< 


2.8 & 10.3 < log(M./M0) 10.5 














-0 


07 ± 0.03 


A/ 


2.4 < C 


< 


2.8 & 10.5 < log(M«/M0) < 10.8 


0.11 ±0.04 


A/ 




0.14 ±0 


08 


A/ 


-0 


09 ± 0.02 


A/ 


2.4 < C 


< 


2.8 & 10.8 < log(M*/M0) < 11 














-0 


13 ±0.02 


A/ 


2.4 < C 


< 


2.8 & log(M./M0) > 11 


-0.12 ±0.04 


A/ 










-0 


21 ±0.05 


A/ 


C ^ 2.4 


& 


10.3 < log(M*/M0) ^ 10.5 








0.23 ±0 


15 


A/ 








C < 2.4 


& 


10.5 < log(M./M0) 10.8 




















C < 2.4 


& 


10.8 < log(M./M0) < 11 






-0.39 ± 0.11 A/ 


-0.17 ±0 


05 


A/ 


-0 


12 ± 0.06 


A/ 


C ^ 2.4 


& 


log(M./M0) > 11 








-0.07 ±0 


12 


A/ 


-0 


04 ± 0.02 


A/ 



Table 3. Systematic uncertainties on the total stellar mass density (p*), the mass density of metals in stars (pz) and the average stellar 
metallicity {{Zt /Zq)). We first estimated these quantities using the stellar metallicities and stellar masses corrected according to Tables[T] 
and [2] The systematic uncertainties are then expressed as the difference between the 'corrected' and default values. 



Ap.(xlO») Apz(xl06) A{Z,/Zq) 
(1) (2) (3) (4) 



stacking 0.569 : 0.063 2.184 : -0.272 0.125 : -0.060 

a/Fe 0.176 : 0.020 0.815 : -0.437 0.065 : -0.070 

prior -0.249 : -0.488 0.403 : -0.873 0.145 : 0.025 

aperture -0.554 -1.943 -0.14 

total magnitude -0.419 -0.922 -0.010 



There are relatively fewer determinations of the mass 
density of metals in stars at the pre sent epoch. We fin d 
agreement within la with the census bylDunne et al] l|2003h . 



and with the estimate bvlPagell (|2002|) based on the stellar 
■ I]. 



mass density of iFukugita et al.l ( 199S ) and assuming solar 
metallicity for stars. In contrast. iFukugita fc Peebles! (|2004l ) 
derive a metal mass density of only 4.6 x 10^/i70 M0 Mpc~^ 
in main-sequence stars and substellar objects (their table 3), 
i.e. roughly 60 percent of the total metal mass density in 
stars that we obtain. 



Cosmic 



hydrodynamic 



simulations 



of 



iDave fc Qppenheimerl (|2007t ) predict a total metal mass 
density in stars of only 3.4 x 1O^M0 Mpc""^, which is at 
the lower end of the observational range. On the other side, 
integration of dust-corrected SFH data points overpredicts 
the amount of metals presently locked up into stars (as 
it overpredicts the stellar mass density) , if we assume 
pz = yp* with y = 0.023 (|Madau et all 1 19961 ). Assuming 
a lower yiel d of y — 0.015, following a recent estimate by 
IConti et al.l l|2003t ). brings the integration of the cosmic 
chemical enrichment history in better agreement with ob- 
servations. The integration of the cosmic SFH predicted by 
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the models of lPei et al.l ^99^ provides a better agreement 
with observations, but a higher metal yield seems to be 
favoured in this case. A similarly good agreem ent is found 
with the predictions from iCalura fc Matteu ccii (|2004l '). who 
study the mean metal abundances for galaxies of different 
morphological types by means of chemo-photometric models 
for ellipticals, spirals and irregular galaxiesl3 These models 
differ fr om cosmic chemical evolution models (such as 
those bv lPei fc Faiilll995l : iPei et al.|[l999l ') in that they are 
designed to study different galaxy types and abundances 
of single elements, rather than average properties of the 
Universe. Finally, we also calculate the metal density in 
stars predicted by the Millennium Simulation, obtaining a 
value of 4.053 x lO^M© Mpc"^. This is ~ 1.5a lower than 
our SDSS-based result (see Section [4.311 . 

As a general remark, it seems that hydrodynamical sim- 
ulations have greater difficulties (with respect to chemical 
evolution models) in predicting the total metal abundance, 
probably due to the sensitivity to the feedback efficiency. 
Indeed, as far as the averag e stellar metallicity is concerned, 
while the compilation from lCalura fc Matteuccl (|2004l ) pro- 
vides a mean metal abundance in stars in agreement with 
our estimate, recent hydrody namical simulati o ns see m to 
predict lower values of (Z*). iKobavashi et all (|2006 l) pro- 
pose hypernova feedback as an extra source of feedback to 
reach agreement between their simulations of cosmic chemi- 
cal enrichment and the observed fraction of baryons in stars 
(less than 10 percent; we find a stellar baryon fraction of 7 
percent). However, this prescription would provide average 
stellar metallicity of about 0.7 Zq, lower than our value, 
even accounti ng for systematic uncert ainties. Similarly, sim- 
ulations from IPave fc Qppenheimerl (.2007 ). which include 
a self-consistent treatment of enriched outflows, predict a 
mean stellar metallicity today of half solar, roughly 3a lower 
than what we derive. In their simulations, the star-forming 
gas in galaxies reaches instead a mean metallicity slightly 
above solar by the present epoch. A higher average metal- 
licity in the star-forming gas with respect to stars might be 
expected because the gas-phase abundance traces the abun- 
dance of the last (most metal-rich) generations of stars and 
not the average abundance of the entire stellar population. 
Com bining gas-phas e oxyg en abundances of SDSS galaxies 
from Tremonti et al.l (|2004l ) , weighting each galaxy by its gas 
masj^'^l and l/Vmav. we find an average gas-phase metallic- 
ity of 12 -I- \og{0/H) — 8.93. However, we cannot robustly 
assess that the gas-phase metallicity is indeed higher than 
the stellar one, b ecause of uncertainti es in the solar oxygen 
abundance scale (|Asplund et al.ll2005l ). and a potential sys- 
tematic overest imate introdu c ed b y strong-line metallicity 
indicators (e.g. iBresolin et al.ll2004i' ). 



Table 4. Present-day mass density of metals in stars derived 
from this work and from the literature. All values have been 
adjusted to the concordance A CDM cosmology with Ho = 
70 h7o km s-i Mpc"!, with p^rit = 1-36 X 10"/ifo Mq Mpc'^, 
and to a Chabrier IMF when necessary. 



Stellar metallicity density 


pz{106 h70 Mq Mpc-3) 


f7z{10-5/l-l) 


Reference 


7 nqq+^ l*^ 

8.96 
4.62 ± 1.08 


6.58 
3.4 ±0.8 


This work 
1 
2 
3 


6.57 
3.4 
4.2 - 6.4 
8.7- 13.3 


4.83 
2.5 

3.09 - 4.7 
6.3 - 9.7 


4 
5 
6 
7 



O bservat i onal e stim ates of metal mass d ensi ty in stars are from: 
(llPagell l|2002l ). ('2 tounne et al.l l|2003h . ('3 lFukugita fc Peebles! 
ll2004f) . We report al so values of pz predicted from chemo- 
photometric models of lCalura fc Ma tteuccii ||2Q04^ (4), from cos - 
mic hydrodynamic simulations of [ Pavi'fc O ppen heime j (l2007t) 
(5), and by integrating the SFH of lPei et a l. (1999) cosmic chem- 
ical evolution models (6) and the analytic fit of i Cole et al.l fcOOll) 
to the dust-correcte d cosmic star formation history (p*) data of 
ISteidel et al.l lll999h (7), assuming pz = yp* (in both cases the 
two values correspond to y = 0.015 and y = 0.023 respectively). 
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To convert their metals density from Salpeter IMF to Chabrier 
IMF we adopt a factor of 1.05 (instead of 1.75 used for the stel- 
lar mass density), which accounts for the higher metallicities ex- 
pected with a Chabrier IMF. While Salpeter-based mass-to-light 
ratios are higher on average by 1.75 than Chabrier-based ones, 
the metallicity is a factor 0.6 lower due to the fewer massive stars 
(> 8Mq) with respect to Chabrier. 

The g as masses are c o mput ed from the SFR surface density 
following pTremonti etal] ||2004|) . 



Figure 5. The estimates of metal and baryon density in stars 
(PZi P*) derived in this work (solid line and grey-shaded region) 
are compared to other observational determinations in the litera- 
ture (blue symbols) and to model predictions (red symbols). The 
number on the x-axis gives the reference according to the caption 
of Tables Hand El 
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Table 5. Present-day stellar mass density as derived in this work 
and in the literature. All values adjusted as described in Table|4] 



Stellar 


mass density 




p. (10** hro Mq Mpc-3) 


n,(io-^ h-^) 


Reference 


Q A^ Q+0.569 
3.21 ±0.47 


Z.OU»_Q 

2.36 ± 0.35 


This work 
1 


3.78 ± 0.44 
1.94 - 4.27 


2.78 ±0.32 
3.091^-22 
1.43 - 3.14 


2 
3 
4 


3.14 ±0.94 

^•'"-0.19 
2.22 ±0.2 
r, Q7+0.74 

3.50 ±0.17 


2.3 ± 0.7 
1.29 ± 0.14 
1.63 ±0.15 

9 11+0.54 
^•J^ 1^-0. 26 

2.57 ±0.13 


5 
g 

7 

8 
9 


4.30 ±0.56 


3.16 ± 0.41 


10 


4.89 ±0.95 


3.6 ± 0.7 


11 


2.80 
2.05 
2.38 

2.86 - 3.81 
1.90 - 5.81 


2.06 
1.51 
1.75 
2.1 - 2.8 
1.40 - 4.27 


12 
13 
14 
15 
16 



Observational es tima tes of stellar mass de nsity are from: (1) 
[Coi^rai] llooj), f2') iKochanek et al.l 1I2OOIII . (3') iFukugita et al.| 
||199S ) from the compila t ion oflPage l' hOOt). (A) 'Glazebrook et al.| 
]2003l'). C5) iBell et al.l ll2003l'). (&i Rudnick et al. (2003), (7) 
iRudnick et al.l ||2006^ fS) [Panter et al.l (2006), (' 9) iDriyer et al.l 
l|2007a^ . do) iDriver et al.l l|2007bl ). (11) Shankar et al.l (12006 



Calura . 



from (12) 

Matteuccil ll2004h . (14) 



Model predictions are 



iPei et al.l 



^ , ^ ^. (13) 

iMonaco et aj] (|2003), (15) 
NagamineetalJ (12006 ) and i ntegration of the a nalytic fit of 
Cole et al.l l|2001^ to p* data of lSteidel et al.l (|l999l ) (16, the two 



values are without and with correction for dust extinction, respec- 
tively) . 



4 THE DISTRIBUTION OF METALS AND 
BARYONS IN THE LOCAL UNIVERSE 

We now study how metals and baryons are distributed ac- 
cording to different galaxy properties (Section l4.1|l and com- 
ment on the characteristic age of the mass and metallic- 
ity distributions as a function of stellar mass (Section 14. 2|) . 
Finally, we compare the observed distributions in detail 
with those predicted by the galaxy formation model of 
iDe Lucia fc Blaizoj (12007 ) applied to the Millennium Sim- 
ulation (jSpringel et al.ll2005l . Section Ojl . 



4.1 



An inventory of the stellar metallicity and 
stellar mass 



In addition to quantifying the total metal budget in stars 
of the local Universe, it is of interest to investigate the 
fractional contribution to the total amount of metals and 
baryons in stars today by galaxies with different properties. 
Which are the galaxies that contain the bulk of the metals, 
and how do they differ from the galaxies that contain the 
bulk of the stellar mass in the local Universe? In order to 
answer these questions we plot in Figs. [6] and [7] the fraction 



of the total mass of metals in star^3 a function of various 
galaxy properties. In Fig.|6]we analyse observable quantities, 
such as the concentration parameter (a), the rest-frame g — r 
colour (b), the 4000 A- break index strength (c), the stellar 
velocity dispersion (d) and the absolute r-band magnitude 
(e). The distribution as a function of the derived physical 
properties is shown in Fig. [T] r-band light-weighted age (a), 
mass- weighted age (b), stellar mass (c) and stellar metallic- 
ity (d) . The grey-shaded histogram gives the distribution for 
the sample as a whole, while the dotted and the dot-dashed 
lines represent the contribution from the high-S/N galaxies 
only and from the low-S/N galaxies only (as derived from 
the stacked spectra), respectively. It is evident that, neglect- 
ing low-S/N galaxies, we would have missed a substantial 
fraction of the total amount of metals in the local Universe, 
in particular at low velocity dispersions, low concentrations, 
low D4000 values and hence young ages, while there is no 
strong segregation in luminosity and stellar mass. 

The red solid line in each panel traces for comparison 
the fraction of the total stellar mass as a function of the 
different parameters. The stellar mass density distribution 
for SDSS galaxies has been studied as a function of spec- 
tral and photometric properties of galaxies, of their size 
and morphology, stellar mas s and surface mass dens ity by 
iKauffmann et al] l|2003l ) and iBrinchmann" et al.l(|2004l ). Our 
distributions agree with those previously derived, although 
some differences may be expected due to the different sample 
definition. In particular, notice that the stellar mass distri- 
bution as a function of the concentration parameter (Fig. [6^) 
is strongly double-peaked, where as the distribution shown 
by IKauffmann et al.l (|2003l ) and IBrinchmann et all (|2004l ) 
does not peak at any particular value. This is likely an arti- 
fact in our distribution caused by the definition of average 
concentration for the coadded spectra (see Fig. [T};). In addi- 
tion to the parameters already studied we are able to show 
here the distribution of stellar mass directly as a function of 
age and not only of D4000. 

Thanks to the good statistics provided by the SDSS 
DR2 we can give accurate description of the distributions 
shown in Figs. |S] and [T] For the velocity dispersion log ov 
and the absolute r-band magnitude we are limited by the 
size of the bins in which low-S/N galaxies are grouped to 
obtain high-S/N coadded spectra (0.05 dex and 0.5 mag 
respectively). In Table |6] we give the mode of each distri- 
bution, which indicates the typical parameter of the galax- 
ies where metals are most likely found (last column). The 
5, 10, 25, 50, 75, 90, 95 percentiles of each distribution are 
listed as well. The same quantities are given also for the dis- 
tribution in stellar mass density. In the last row of Table (HI 
we indicate the fraction of the total stellar mass contained 
in galaxies that contribute different fractions of the total 
metal content. The systematic uncertainties have been es- 
timated by calculating the distributions with the masses, 
metallicities and ages corrected following Tables [1] and (2] 
The uncertainties quoted here give the range of variation in 
the distributions. 

From Figs. |6] and [7] and Table |6] it appears that the 
distribution of the metals locked up in stars does not differ 



1'* The differential pz normalized by the total density of the 
whole galaxy sample. 
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substantially from the distribution of the stellar mass. 
In other words, the galaxies that contribute a significant 
fraction of the total amount of metals in stars today are 
also those that contain most of the total stellar mass. 
The similarity in the two distributions is determined 
by the relatively narrow range (roughly two orders of 
magnitude) in mass covered by the stellar mass distribu- 
tion: most of the weight is concentrated in galaxies with 
stellar mass around IO^^Mq. The particular shape of the 
mass-metallicity relation does not influence substantially 
the stellar metallicity distribution. However, the increase 
of metallicity with mass causes the stellar metallicity 
distribution to be shifted to slightly higher values of stellar 
mass (see Fig. [T];). Indeed, differences are more clearly 
evident in galaxies with low concentrations, low velocity 
dispersions and young ages: they contain a non-negligible 
fraction of the total stellar mass (almost comparable to 
that contained in high-concent ration galaxies), becaus e 
they are more numerous (see e.g. iBrinchmann et al.ll2004 ). 
but their stars contain a much smaller fraction of met- 
als. We describe these results in more detail in the following. 

We can characterize the properties of the typical 
galaxy contributing stellar mass or metals by the median 
of the corresponding distributions. It is remarkable that 
the typical galaxy in the local Universe appears to have 
global properties close to that of the Milky Way and 
M31. More quantitatively, it has a velocity dispersion 
of ~ 130 km s~^, an absolute magnitude about 1 mag 
brighter than M* = —20.28 in r-bancF^. a, g ~ r colour in 
agreement with what expected from the colour-magnitude 
relation of elliptical galaxies, a concentration parameter of 
2.7 (characteristic of a galactic disk with a significant bulge 
component), D4000 of 1.7 corresponding to a fairly old 
stellar age of ~6 Gyr, and a typical mass of ~ 6 X 10^°Mq. 

Fig. [6^ shows the distribution of mass and metals as 
a function of the concentration parameter. The peak of 
the distribution of metals is at a concentration parameter 
of ~2.9 (therefore it is contributed by galaxies that are 
predominantly early-type). The fraction of metals drops 
quickly in galaxies with concentration parameter below the 
median value of ~2.7. While early-type galaxies (C ^ 2.8) 
contain roughly 40 percent of the total metal budget 
in stars, late-type galaxies (C ^5 2.4) contribute less 
than 25 percent. The contribution of early- and late- type 
galaxies to the total metal budget in stars increases to 
60 and 40 percent respectively if C = 2.6 is adopted 
as thresh old to separate lat e- and early-type galaxies 
(following Istrateva et al.l 1200 ll ) . We note that the c hemo - 
spectrophotometric models of ICalura fc Matteuccil (|2004l ) 
predict that, while spheroids are the largest contributors 
to the total amount of metals (in different phases) in the 
present Universe, they also contribute significantly to the 
enrichment of the IGM and the majority of metals in stars 
come instead from spiral galaxies (60 percent against the 
40 percent contributed by spheroids). As far as the stellar 
baryon fraction is concerned, we find that the two classes 
of galaxies contribute the same fraction of the stellar mass 

15 From lBlanton et al.l l l2003h . corrected to z = 0. 



density (30 or 50 percent, depending on which of the two 
concentration cuts we adopt). This is consiste nt with what 
was already found bv lKauffmann et al.l ()2003l ). 

The distributions of metals and baryons are very 
similar to each other also as a function of colour (Fig. [BJa). 
Red galaxies contribute the same fraction to metals as 
to baryons in stars. The metal fraction becomes smaller 
than the stellar mass fraction only in galaxies bluer than 
g ~ r — 0.5. At least half and up to 75 percent of the total 
stellar mass (and metals) are contained in red-sequence 
galaxies (assuming g — r = 0.7 or 0.6 as colour cut respec- 
tively). This result is in good agreement with what was 
found bv lBell et al.l (|2003l ) separating elliptical galaxies with 
a ma gnitude-dependent colour cut, and by iBaldrv et al.l 
faooi) who also distinguish between red-peak galaxies 
according to their distribution in the colour-magnitude 
plane. We note that these fractions are also consistent with 
the distributions as a function of concentration discussed 
above, given that 84 percent of C ^ 2.6 galax ies satisfy also 
the colour-based selection of iBell et al.l (|2003l ). 

The distribution of stellar metallicity as a function of 
D4000 (Fig. [BJ:) deviates from the distribution of stellar 
mass in the range of D4000 occupied by late-type, star- 
forming galaxies. Both distributions show a strong peak at 
D4000=1.9: galaxies with D4000 above this value contain 
roughly 25 percent of the metals and 25 percent of the mass 
in stars today. The distribution in mass is clearly bimodal 
and has a secondary peak at D4000 ~ 1.4. Galaxies with 
D4000 1.4 contribute another 25 percent to the total 
stellar mass, but only 10 percent of the metals. 

Fig. [7] complements the picture derived from obser- 
vational properties with physical ones. Panels (a) and (b) 
illustrate that the differences in the stellar metallicity and 
stellar mass distributions with respect to D4000 and colour 
are reflected in the stellar age. Focusing on the r-band 
light- weighted age (Fig. [7^), galaxies older than 8.5 Gyr 
contribute the same fraction (25 percent) of the total stellar 
mass and the total stellar metallicity densities in the local 
Universe, and only 5 percent comes from galaxies older 
than 10 Gyr. The distribution in stellar metallicity declines 
rapidly at ages younger than 6.3 Gyr {log{tr/yr) ~ 9.8), 
where roughly 50 percent of the total stellar mass, but less 
than 40 percent of the total amount of metals, comes from. 
Similar results are found considering the mass-weighted 
age (Fig. EJd). The only significant difference is that the 
distributions here are narrower due to the older mass- 
weighted age in young, low-mass galaxies with respect to 
their light-weighted age. 

The dependence of the fraction of mass and metals 
in stars on stellar mass is shown in Fig. [7J;, quantifying 
what expected on the basis of the distributions against 
velocity dispersion and absolute magnitude, both tracers 
of the total stellar mass, especially in quiescent elliptical 
galaxies (Fig. [6ji,e). Half of all the metals locked up in 
stars today are contained in galaxies more massive than 
7.2 X IO^^Mq (or with velocity dispersion higher than 
148 km s~i), which contain roughly 40 percent of the 
total stellar mass. Galaxies with masses below 4 X 10^°Mo 
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(or av < 110 km s~^) contain only 25 percent of the 
total metal budget and about 35 percent of the total 
stellar mass. Similarly, panel (d) shows the (mass- and 
volume-weighted) projection of the mass-metallicity rela- 
tion onto the metallicity axis. This illustrates, consistently, 
that at least half of the metals are contained in galaxies 
with metallicity above solar, which are predominantly 
massive ellipticals and the bulges of massive late-types, 
galaxies with masses above IO'^^Mq. The steepening of 
the mass-metalhcity relation becomes clear at metallicities 
below 0.7 X Zq (or below the transition mass of lO^^ '^Mo), 
where the largest differences in the relative contribution to 
the amount of baryons and to the amount of metals are seen. 

In conclusion, we find that the bulk of the total metals 
locked up in stars in the local Universe resides in galaxies 
with masses just above the transition mass in the mass- 
metallicity relation, with morphology and spectral proper- 
ties of intermediate-type galaxies (early late-types or ellipti- 
cals), and with fairly old stellar populations. Given the shape 
of the mass density distribution, these results are in agree- 
ment with the correlations between stellar metallicity, age 
and stellar mass studied in paper I. Nonetheless, late-type, 
star-forming galaxies (with masses below the characteristic 
mass of 3 X 10^" Mq, low D4000 values and concentration 
parameters characteristic of disc-dominated galaxies) con- 
tribute roughly 20 percent of the total mass density of met- 
als in stars and a slightly higher fraction of the total stellar 
mass density (30 — 35 percent). 

4.2 The characteristic age of the mass and 
metallicity distributions 

We investigate further the distribution of metals and 
baryons as a function of stellar age. As shown in Fig. [7^,b 
both the steUar mass distribution and the metals distribu- 
tion have a peak at a mean age of \og{t/yr) — 9.95 (al- 
most 9 Gyr). If we could translate this characteristic age 
into a redshift, it would correspond to a characteristic for- 
mation redshift of z ~ 1.4, which nicely falls in the red- 
shift range over which the cosmic metal production rate and 



Lillv et al.lll996: Madau ct al. 


1996. 1998:lGlazebrook et al.1 


20041: HoDkins & Beacomll200a 


, and references therein) . This 



is of course only a naive interpretation, because we can only 
assign an average age (with a larger weight to younger stars) 
to individual galaxies and we cannot describe the real dis- 
tribution of stellar ages within a galaxy. This is even more 
critical when the distribution in stellar metallicity is consid- 
ered, as long as we assign a fixed metallicity to all stars in 
a galaxy, rather than following the chemical evolution along 
the star formation history. It is however interesting to notice 
that the result obtained from this 'archaeological' approach 
is close to what expected from 'direct' investigation of the 
cosmic star formation history. 

Moreover, there is evidence that the timescale of star 
formation depends on the mass of the galaxy, with more 
massive galaxies having an early and shorter star forma- 
tion, and less massive galaxies having a star formation more 
extended toward the present day. If so, the distribution of 
stellar ages in individual massive galaxies would be narrower 
than in less massive galaxies around the (mass- or light- 



weighted) average age, which would hence be more repre- 
sentative of the average formation redshift of the stars. It is 
thus useful to look at the distributions of baryon and metal 
densities in stars as a function of stellar age for galaxies with 
similar masses. These are shown in Fig. |8] (left- and right- 
hand panel respectively). The black continuous line shows 
the distribution of p* and pz for the sample as a whole, while 
different colours and line styles correspond to different stel- 
lar mass bins (each distribution is normalized to the total 
p, and pz in the corresponding mass bin). For reference, we 
indicate in the upper x-axis the redshift corresponding to 
the age on the lower x-axis (interpreted as lookback time 
from the present). 

The distribution for the entire galaxy population is 
traced very well by the distribution of galaxies in the mass 
range between 3 x 10^" and IO^^Mq, which alone contain 
38 percent of the total stellar mass budget (these galaxies 
also provide the largest contribution to the star formation 
rate d ensity and ste l lar m a ss density up to z ^ 1, as shown 
e.g. bvlPanter et al.1 l|2006l ): lBorch et al.l l|2006h : IZheng et all 
(j2007|)). This is true for both the stellar mass density and 
the metal density in stars. The contribution to the total 
baryon and metal budget in stars from galaxies with old 
stellar populations increases significantly from the lowest- 
to the highest-mass bin. This refiects the median relation 
between the galaxy mean age and stellar mass (see e.g. fig. 8 
of paper I). Not only do the distributions in mass- weighted 
age shift gradually to younger ages as less massive galaxies 
are considered, but they also span larger ranges in agj^. a 
result of the increasingscatter in age at lower masses in the 
age-mass relationship!^ 

Quantitatively, 25 percent of the total stellar mass and 
metal density of the most massive galaxies (> lO^^Mo) 
is contained in galaxies older than 9.5 Gyr, and half of it 
in galaxies older than 8.5 Gyr (corresponding to redshift 
greater than 1.2). As much as 75 percent of the stellar mass 
and metal density of low-mass galaxies (below 10^" Mq) 
is distributed in galaxies with mass-weighted ages younger 
than ~5.5 Gyr, and 50 percent in galaxies younger than ~4.5 
Gyr (corresponding to redshifts below ~0.5). This trend is 
illustrated in Fig. [9] In each of the stellar mass bins de- 
fined in Fig. [S] we define two characteristic mass- weighted 
ages as the median of the distribution of stellar mass and 
the distribution of metals versus age. These characteristic 
ages are plotted as a function of stellar mass in Fig. (5] We 
also indicate the effect of individual sources of systematic 
uncertainties in each mass bin. The aperture bias and the 
scaled-solar abundance ratio of the models, in particular, 
affect to a larger extent more massive galaxies. 

Despite the large uncertainties, it is clear from Fig. [9] 
that the characteristic ages of the stellar mass and metallic- 
ity distributions become progressively younger in less mas- 
sive galaxies. This is in agreement (at least qualitatively) 



It is interesting to notice that the bimodal distribution ex- 
pected from the bimodality in D4000 appears in galaxies of inter- 
mediate masses. This would be even clearer if we plotted light- 
weighted age instead of mass-weighted age. The bimodality in age 
is instead almost lost when the full population is considered. 

Note that these are logarithmic ranges, which means that it is 
the relative age range, not the absolute one, which increases on a 
linear scale toward lower masses. 
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Figure 6. The fraction of the total mass of metals locked up in stars in the local Universe is shown as a function of various observable 
galaxy properties: (a) the concentration parameter, (b) the rest-frame g — r colour, (c) the 4000A-break index strength, (d) the stellar 
velocity dispersion and (e) the absolute r-band magnitude. The dotted line shows the contribution from high-S/N galaxies, while the dot- 
dashed line shows the contribution from low-S/N galaxies (obtained from high-S/N stacked spectra, see text for details). The continuous 
black line (grey-shaded histogram) shows the distribution obtained when both contributions are taken into account. The distribution 
of stellar mass as a function of the various properties is described by the red solid line in each panel. Note that the resolution in the 
distribution versus velocity dispersion and absolute magnitude is limited by the width of the bins in which low-S /N galaxies are grouped 
to obtain the co-added spectra (0.05 dex and 0.5 mag respectively). 



with the findings that the characteristic epoch at which the 
bulk of the stars (and metals) in present-day massive (mostly 
elliptical) galaxies formed is at redshift at least greater than 
1 (e.g. Irhomas ct al. 2005, and references therein). Finally, 
it is worth noting that the characteristic age of the distri- 
bution of metals is never younger than the characteristic 
age of the distribution of mass. This is particularly true 
at masses below 10^^ M0. This probably indicates that the 
correlation between age and stellar metallicity persists in 
individual stellar mass bins, with older galaxies being also 
more metal-rich. This correlation is weaker in higher-mass 
galaxies, because of the narrower ranges in both age and 
metallicity and because contrasted by the age-metallicity 
degeneracy (see e.g. figs. 11 and 12 of paper I). 

4.3 Comparison with the Millennium Simulation 

The distribution functions of mass density of baryons and 
metals in stars provide important quantitative constraints at 
redshift zero against which models of galaxy formation and 
evolution in a cosmological context should be tested. We 



provide here a first comparison between our observational 
results and model predictions. We consider the results of the 
Millenmurn Run, the largest N-body simulation of structure 
formation carried out so far within the ACDM cosmological 
paradigm HSpringel et al.ll2005l ). We make use of the galaxy 
catalogues produced by semi-analytic galaxy formation al- 
go rithms run on the M illen nium merger trees, as desc ribed 
in lCroton (|2006l ) and lDe Lucia fc Blaizo3 (|2007l ). and 

select galaxies at 2 = 0. 

Before going into the detail of the distributions of met- 
als and baryons in stars, we calculate the total p* and pz 
predicted by the simulation. The total stellar mass density 
is 3.203 X 10* M© Mpc~'^, in agreement with our determi- 
nation as already mentioned in Section [3.31 The total mass 
density of metals in stars is instead 4.053 x 10^ Mq Mpc"'', 
lower than what observed (in particular it is 1.5(j lower than 
our result). This originates from the narrow metallicity dis- 
tribution predicted by the models and the lower metallicity 
for massive galaxies, as we discuss below. 

It is also of interest to mention the metal and baryon 
mass densities in other gaseous components, inside or out- 
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Figure 7. Same as Fig. |6] but here the fraction of stellar mass and the fraction of metals in stars are shown as a function of the derived 
physical parameters: (a) r-band luminosity weighted age, (b) mass- weighted age (both ages are corrected to 2 = by the lookback time), 
(c) stellar mass and (d) stellar metallicity. 




Figure 8. Differential distribution of stellar mass density (left panel) and mass density of metals in stars (right panels) as a function of 
mass- weighted age for the whole sample (continuous line). Different line styles distinguish the distribution for galaxies in different mass 
bins. Each distribution is normalized to the total p« and pz in the corresponding mass bin. For reference the upper x-axis gives the 
redshift corresponding to a given mean age. 
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Table 6. For each parameter X, the percentiles of the distribution of stellar mass and of metals as a function of X are given. In the 
last column we indicate the mode of each distribution. The last row quotes the fraction of the total stellar mass in the local Universe 
contained in galaxies that contribute different fractions of the metal content. 
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fThe estimated systematic uncertainty is smaller than 0.01, i.e. less than 20 percent of the bin size. 



side galaxies, as predicted by the Millennium Simulation. 
The mass and metal densities in stars represent the 17.6 
and 25.3 percent of the corresponding total values. Sim- 
ilar contributions (13.5 and 11.8 percent for the baryons 
and the metals, respectively) come from the gas ejected 
fro m galaxies into the s urrounding intergalactic medium 
(see lPe Lucia et al.ll2004h . The cold gas inside galaxies pro- 
vides only 6.6 percent and 9.8 percent to the total den- 
sity of baryons and metals, respectively. The remaining 62.2 
and 53.1 percent is contained in the hot gas component. 
The baryon fractions derived from the Millennium Simula- 
tion are in agreement with those predicted by th e chemo- 
photometric models of ICalura fc Matteuccil (12004! ) ■ In both 
cases the fraction of baryons i n stars is higher than what pre - 
dicted by the simulations of iDave fc Oppenheimerl (|200'7| ) 
and than what we observe. The mass fraction of metals 
in stars predicted by the Mille nnium Simulation is i nstead 
abo ut half of wha t quoted by ' Calura fc Matteuccil (|2004l ) 
and lDave fc Qppe nhcimcr (2007). 

We now discuss the contribution to the mass and metals 
in stars by different galaxies. After selecting galaxies at z — 
from the Millennium Catalogue, we construct the stellar 
mass-weighted distributions of baryons and metals in stars 
as functions of the galaxy stellar mass, mass- weighted age 



and stellar metallicity. Fig. [TD] compares such distributions 
(gray-shaded histograms) with those derived from our SDSS 
measurements (solid lines). 

Panels (a) and (d) of Fig. [TD] illustrate the fraction of 
mass and metals as a function of stellar mass. The distri- 
butions predicted by the models display an offset in stellar 
mass with respect to the observationally derived distribu- 
tions of ~0.06 dex. We note that this offset is reduced when 
we use stellar masses based on Petrosian instead of model 
magnitudes (as shown by the dotted line). Comparison be- 
tween the shaded and dotted histograms in panels (a) and 
(d) shows a good agreement between the observational re- 
sults and the predictions from the simulations. It is clear, 
however, that the models predict too much mass in low- 
mass galaxies with respect to observations. This feature of 
the models has been already observed in the comparison 
with observational data of the bj- and K-ba nd luminosity 
funct i on and the stellar m ass function (see e.g. Croton et al.l 
I2OO6I : iBertone et ai]|2007^ . Recently, iBertone et al.1 (|2007D 
have presented a new implementation of a feedback scheme 
which incorporates a dynamical treatment of galactic winds 
powered by supernovae and by stellar winds, as alternative 
to empirically-motivated schemes of supernovae feedback, 
but without replacing AGN feedback. This treatment seems 
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Figure 9. Mass-weighted age above which half of the total 
amount of mass (empty circles) and metals (filled circles) in stars 
are contained at the present epoch, as a function of stellar mass. 
Galaxies have been divided in the same five bins of stellar mass 
as in Fig. [8] Boxes and arrows quantify the ranges of the different 
sources of systematic uncertainties. 

to alleviate the discrepancy in the galaxy stellar mass func- 
tion, by reducing the abundance of low-mass galaxies and 
reaching a better agreement, in this respect, with observa- 
tions. 

Panels (b) and (e) show the fraction of baryons and 
metals in stars as a function of the galaxy mass-weighted 
age. While the predicted and observed distributions agree 
reasonably well at the oldest ages, there is clearly a deficit 
of young galaxies (with ages younger than logt = 9.8) in 
the simulations. This discrepancy could originate, at least 
in part, from the fact that, although we use mass-weighted 
ages for this comparison, the observationally derived mean 
ages are always extracted from fits to the galaxy spectra, and 
hence are always somehow weighted by light. This could thus 
give more weight to the youngest stellar populations in the 
galaxy. We note however that the tendency of the models 
to produce too old stellar populations with respect to those 
that we derive appears also in the relation between the light- 
weighted age and the stellar mass o f elliptical galax i es (as 
seen in the comparison of fig. 6a of lDe Lucia et al.l (120061 ) 
and fig. 17d of paper II), even though there is a good qual- 
itative agreement. 

Finally, panels (c) and (f) compares the predicted and 
observed distribution as a function of stellar metallicity. The 
disagreement between the observational result and model 
prediction is caused by the narrow stellar metallicity range 
spanned by the simulations. Moreover, while the observa- 
tional trend of increasing metallicity with increasing mass 
is reproduced by the models, the predicted relation is too 
shallow with respect to observations, hardly reaching so- 
lar metallicity for the most massive galaxies. The galactic 



wind feedback scheme implemented bv lBertone et al.l (|2007l ) 
seems to provide a better agreement with observations, caus- 
ing a steepening of the mass-metallicity relation (due to 
the high efficiency of mass ejection in small haloes which 
suppresses star formation and hence the amount of metals 
that can be produced and locked up in stars) and predicting 
super-solar metallicities in massive galaxies. 

As a general remark it is also interesting to notice the 
similarity between the distributions of mass and the cor- 
responding distributions of metals in stars. While this is 
true also for the observationally derived distributions (as 
discussed in Section 14. ip , it appears to a greater extent in 
the models. This originates from the narrower range in stel- 
lar metallicity and the shallower stellar metallicity/stellar 
mass relation with respect to observations. 



5 SUMMARY AND CONCLUSIONS 

In this work we have exploited recent estimates of physical 
parameters, such as stellar metallicity and stellar mass, for 
a comprehensive sample of more than 10^ nearby galaxies to 
derive the total mass density of metals and baryons locked 
up in stars in the local Universe, also expressed in terms 
of the present-day mass- weighted average stellar metallicity. 
Moreover, we have quantified the contribution to the total 
amount of metals by galaxies with different morphological 
and spectral properties, and compared the distribution of 
the stellar metallicity density with that of the stellar mass 
density. 

The sample used is drawn from the SDSS DR2 and it 
includes galaxies spanning a wide range of star formation 
activity, from quiescent early-type to actively star-forming 
galaxies. The stellar metallicities, ages and stellar masses 
of these galaxies were previously derived (paper I) by com- 
paring an optimally selected set of spectral absorption fea- 
tures to a comprehensive Monte Carlo library of model 
star formation histo r ies (S FH), based on the high- resolution 
iBruzual fc CharlotI l|2003l ) population synthesis code. We 
have shown in paper I that the uncertainties on the derived 
physical parameters directly depend on the mean spectral 
signal-to-noise ratio (S/N). Because of this, in previous work 
we have focused on galaxies with S/N greater than 20, bi- 
asing our sample towards more concentrated, higher surface 
brightness and higher velocity dispersion galaxies. 

We need to include here also galaxies at lower S/N in 
order to derive a fair estimate of the total metallicity of stars 
at z = 0.1. We do this by coadding (weighting by 1/Vmax) 
the spectra of low-S/N galaxies with similar r-band abso- 
lute magnitude, velocity dispersion and D4000 until a min- 
imum S/N requirement is satisfied. This allows us to derive 
estimates of the (1/Vmax-weighted average) stellar metal- 
licity of low-S/N galaxies with an uncertainty not greater 
than 0.2 dex, i.e. comparable to the accuracy with which 
stellar metallicity is derived from individual high-S/N spec- 
tra. A similar improvement is obtained for light- and mass- 
weighted age and stellar mass, although these parameters 
are less affected than stellar metallicity by the quality of the 
spectrum. 

We estimated the total mass density of baryons and 
of metals locked up in stars in the local Universe, by com- 
bining the contribution of the individual high-S/N galax- 
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Figure 10. Comparison between predicted and observed distributions of stellar mass and metals. The fraction of metals in stars (upper 
panels) and the fraction of mass in stars (lower panels) are shown as functions of stellar mass (a,d), mass- weighted age (b,e), stellar 
metallicity (c,f). The shaded histograms show the distributions derived using the semi-analytic galaxy catalogue of the Millennium 
Simulation. The continuous line shows the results from our analysis on SDSS. The dotted histogram in panel (a) and (d) shows the 
results obtained when normalizing the mass-to-light ratio by the ^-band Petrosian magnitude instead of model magnitude to derive 
stellar mass. 



ies and the low-S/N galaxies, as derived from the coadded 
spectra. We find, respectively, p. = 3.413 ± O.OOSlfjjs! x 
10* hro Mq Mpc"^ and pz = 7.099 ± O.OWtlUt x 
10*^ hro Mq Mpc"^ This is in broad agreement with other 
measures from the literature, which however span a rela- 
tively large range. The measure of the metal mass den- 
sity is perhaps more significant in discriminating different 
model predictions or approaches. We find a reasonably good 
agreement with resu l ts of c he mical evolution mode l s, suc h 
as those of |Pei et all l|l999l) or ICalura fc Matteuccl (|2004l ). 
Hydrodynamic simulations seem to predict values of pz at 
the lower end of the observational constraints. The total 
amount of metals locked up in stars is sensitive to the feed- 
back efficiency, which regulates the amount of gas and metals 
expelled from galaxies. Integration of cosmic star formation 
histories, on the other hand, in general overpredicts the to- 
tal amount of metals in stars. This is probably due to poor 
knowledge of dust corrections at high redshift and uncer- 
tainties on the average metal yield to convert star formation 
rate into metal production rate. 

By combining the densities of mass and metals in stars 
we can estimate the average stellar metallicity today to be 
consistent with solar (Equation |6} . This is not surprising, 
given that the stellar mass budget is dominated, at least 
in the local Universe, by galaxies just above the transition 
mass (approximately L* galaxies) , which have typically solar 
metallicity. This value is in agreement with the prediction 



of the chemo-photomet ric models of [Calura fc Matteuccil 
(2004') and the result of I Edmunds fc PhillippsI (|l997f r 

The good statistics available from the SDSS allows us to 
provide an analysis of the distribution of the stellar metal- 
licity density as a function of global galaxy properties (such 
as velocity dispersion, luminosity and stellar mass), as a 
function of morphology (as approximated by the concen- 
tration parameter), and as a function of spectral and physi- 
cal properties of the stellar populations (such as colour and 
4000A-break strength, stellar age and metallicity). We have 
compared such distributions with the corresponding distri- 
butions of stellar mass density. Our study has shown that 
the stellar metallicity and the stellar mass distributions do 
not differ significantly, in particular at the high-mass end. In 
other words, the galaxies that contribute most of the total 
amount of metals in stars have properties similar to those 
containing the bulk of the total mass in stars. 

The typical galaxy in the present-day Universe where 
baryons and metals reside has a concentration parameter 
of 2.7 (indicating a disc and a significant bulge compo- 
nent), a D4000 = 1.7, corresponding to a fairly old stel- 
lar age of ^6 Gyr, and a stellar mass of 6 x IO^^A^q. Note 
that the characteristic mass found here is above the stellar 
mass of 3 X 10^°Mq at which the transition from low-mass, 
metal-poor, disc-dominated galaxies to high-mass, metal- 
rich, bulge-dominated galaxies occurs. Galaxies with masses 
greater than 3 x 10^^ Mq, where the mass-metallicity rela- 
tion starts to flatten (|Tremonti et al.ll2004l : iGallazzi et al.l 
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I2005I ). contain more than 80 percent of the total amount of 
mass and metals in stars. 

The stellar metallicity and stellar mass distributions 
differ most significantly with respect to concentration and 
stellar age (as indicated also by the colour and D4000), in 
the sense that the relative contributions to the total stellar 
mass and to the total metallicity are reversed in late-type 
galaxies, compared to early-type galaxies (morphologically 
or photometrically classified as such). Separating galaxies 
on the basis of their concentration parameter, we find that 
bulge-dominated and disc-dominated galaxies provide the 
same fraction of the total mass in stars (30 or 50 percent 
depending on the cut in C adopted). On the contrary, the 
contribution of disc-dominated galaxies to the mass den- 
sity of metals is lower than that of bulge-dominated galaxies 
(with the same cuts in concentration). For example, defining 
early-type galaxies as those having C ^ 2.8 and late-types 
those having C ^ 2.4, the two classes of galaxies both con- 
tribute 30 percent to the stellar mass budget, but 40 and 
less than 25 percent respectively to the total metal budget 
in stars. 

Similar considerations hold for the distributions as func- 
tions of colour or D4000. We find that at least 50 up to 
75 percent of the stellar mass (and of the metals) is locked 
into red-sequence galaxies (depending on the assumed colour 
cut). These frac tions are in good agree ment with previou s 
results from e.g. iBell et all (|200a ) and .Baldrv et alj l|2004h . 
We show that these observational distributions are physi- 
cally translated into the corresponding distributions in age. 
More than 50 percent of the mass and the metals locked up 
in stars reside in galaxies with ages older than 6.3 Gyr. 

It would be tempting to translate the characteristic age 
of the stellar populations contributing most of the metals 
and of the mass into a characteristic redshift of formation. 
Unfortunately, this cannot be done straightforwardly be- 
cause our stellar ages are averages over all the stellar pop- 
ulations in a galaxy, weighted by mass or by light. Thus 
they are very sensitive to the galaxy SFH and, on average, 
closer to the last significant episode of star formation. Given 
the different timescales of star formation expected in galax- 
ies with different mass, it is reasonable to look at the dis- 
tributions of mass and metals against (mass-weighted) age 
for different stellar mass bins. This shows that most of the 
mass and metals come from galaxies older than a character- 
istic age which becomes progressively younger from massive 
to less massive galaxies. Considering only massive galaxies 
(> 10^^ Mq) at least 50 percent of the mass and metals is 
in systems older than 8.5 Gyr, which would correspond to a 
redshift z > 1.2. This redshift is compatible with the redshift 
range 1 < z < 2 at which the cosmic star formation rate den- 
sity starts to decline (e.g. iHopkins fc Beacomll200^ . This is 
consistent with the evidence that the bulk of star formation 
in massive galaxies has occurred at high redshifts, while low- 
mass galaxies are still in a phase of active star f o rmation (e.g. 
Brinchmann fc EUidlioOol: I Juneau et alll2005l : iBundv et all 



20051 ). although galaxies appear to experience the same rate 



of decline in their global star formation from z ~1 to the 
present independently of mass (Zheng ot al. 200^). 

Comparison of the observed distributions of pz and p, 
with those predicted by the Millennium Simulation shows a 
good agreement in the distributions against stellar mass. 
However, an excess at low stellar masses (roughly below 



3 X IO'^^Mq) is visible, which has been already pointed 
out in the lumino sity or stellar mass functions (see e.g. 
ICroton et akl 120061). The agreement between our observa- 
tions and the models becomes worse when mass-weighted 
age and stellar metallicity are considered. Although success- 
ful in reproducing the general tr end of increasing bot h age 
and metallicity with mass (e.g. iDe Lucia et e O] 120061 ) . the 
models predict narrower ranges in these physical parame- 
ters than observed. 

Our study provides a new determination of the to- 
tal amount of metals locked up in stars today and allows 
us for the first time to derive a quantitative and accu- 
rate description of the properties of the galaxies hosting 
different fractions of the metals and of the stellar mass 
in the local Universe. Such distributions represent impor- 
tant constraints on models of the cosmic star formation 
and chemical enrichment histories. The detailed knowl- 
edge of the distribution of metals, coupled to the stel- 
lar mass distribution, will allow a more direct compari- 
son with predictions from semianalytic models of galaxy 
formation and evolution . Ongoing large reds hift surveys, 
like V V DS (iLe Fevre et a l. 2004), GOOD S ( Dickinson et al.l 
2003al: IVanzella et all |2005| ). DEEP 2 JPavis et all l2o'o3r 



Madgwick et"al] |2003| ). zCOSMOS (|Lillv et al.ll2006l ). will 



make it possible to extend this kind of studies to redshifts 
as high as z ~ 1 and thus, not only to build a more consis- 
tent picture of the cosmic star formation history, but also to 
understand which are the galaxies that most strongly con- 
tributed to its evolution since z = 1. 
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